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TECHNICAL NOTE 3723 


CALCULATIONS CF THE FLOW OVER AN INCLINED FLAT PLATE 

AT FREE-STREAM MACH NUMBER 1 

By Walter G. Vincent!, Cleo B. Wagoner, and 
Newman H. Fisher, Jr. 

SUMMARY 


A numerical solution has "been obtained of the complete equations of 
inviscid compressible' flow for the case of an inclined flat pi ate at free- 
stream Mach number 1. The mixed flow about the lower surface of the plat e 
is found by relaxation solution of a boundary-value problem in the holo- 
graph plane. Considerable preliminary analysis is required by the pres- 
ence of the free- stream singularity, which must be incorporated analyti- 
cally into the numerical work. The methods devised for this part of the 
work may have application in other problems of transonic flow. The super- 
sonic flow on the upper side of the plate is found in the physical pi »tip 
by a standard form, of the method of characteristics. The calculations 
here are carried only as far as the end of the separated region that 
appears on the upper surface near the leading edge. The results, which 
are for an angle of attack of 13 °, show the pressure distribution on the 
lower surface and the detailed flow field about the lower surface and the 
leading edge. 

The results for the flow field show that the large changes of velocity 
that occur near the leading edge are confined to a surprisingly small part 
of the field. The stagnation point on the underside of the plate, for 
example, is found to be only 0.0016 of the chord aft of the leading edge. 
(This is in contrast to a value of 0.05 for incompressible flow about a 
plate at the same angle of attack. ) The average radius of curvature of 
the sonic line as it approaches the leading edge is even smaller. To see 
the details here, in fact, it is necessary to plot the results to a scale 
in which the displacement of the stagnation point aft of the leading edge 
is approximately 50 inches. 

The calculations of the pressure distribution indicate that Guder ley’s 
earlier solution on the basis of the transonic small-disturbance theory 
gives reasonably accurate results even at the present moderately large 
angle of attack. If the small-disturbance values are corrected by a method 
due also to Guderley, the error is almost completely eliminated except at 
the leading edge where the corrected solution still cannot represent the 
true flow. 
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INTRODUCTION 


The problem of the two-dimensional flow over an inclined flat plate 
at free-stream Mach number 1 has been discussed by Guderley in reference 1. 
This problem is of interest as an example of the unsymmetric transonic 
flow that occurs around a sharp-nosed airfoil when the angle of attack is 
large compared with the nose angle. (On a plate, where the nose angle is 
zero, any angle of attack meets this description.) A solution of the 
problem was given by Guderley on the basis of the small-disturbance theory 
for transonic flow. This solution, arrived at by an ingenious analytical 
procedure, gives useful information on the pressure distribution on both 
surfaces of the plate. Owing to the approximations involved in the theory, 
however, the results are limited to angles of attack which, though large 
in comparison with the nose angle, are in themselves still small. More 
important, the details of the flow at the leading edge, which are of 
fundamental interest, do not appear in the analysis. 

The work to be described in the present paper was conceived as part 
of a study of the transonic flow around leading edges. To supply the 
details missing in the small-disturbance analysis, it is proposed here 
to solve the flat-plate problem using the complete equations of inviscid 
compressible flow. The exclusion of viscosity may, it is realized, lead 
to some error at the sharp leading edge. If past experience is any guide, 
however, an inviscid solution will be of value toward understanding the 
actual viscous flow. It will also provide, for one specific case, a 
check of the accuracy of a known small -disturbance solution. 

Since the analytical problem for the complete equations is a formidable 
one, the solution in the present work is carried out largely by numerical 
means. 1 As in Guderley’ s analysis, the mixed flow over the lower surface 
of the plate is found by solution of a boundary-value problem in the hodo- 
graph plane. The only real difficulty here is with the free-stream singu- 
larity, which must be incorporated analytically into the numerical work. 

The numerical procedures themselves are an extension of the work already 
reported in references 3 and 4. The supersonic flow on the upper side of 
the plate is found in the physical plane by the method of characteristics. 

In the present report this part of the calculation is carried only as far 
as the end of the closed region of separation that appears on the upper 
surface adjacent to the leading edge. It is hoped eventually to complete 
the solution to the trailing edge, but owing to circumstances outside the 
work this will not be possible in the immediate future. 

The calculated results, which are for an angle of attack of 13°, show 
the pressure distribution on the lower surface and the flow field about the 

1 The analytical solution of the two-dimensional airfoil problem 
for the complete equations has been treated by Frankl in reference 2, 
which came to the authors ’ attention after the present work was under way. 
Fr ankl gives a general solution in the form of an infinite series of 
Chaplygin solutions. The application of this general solution to a spe- 
clfic boundary-value problem, however, is not discussed. 
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lower surface and the leading edge. In view of the known properties of 
flows near the sonic flight speed (see refs. 5 and 6), the results should 
he applicable for a range of free-stream Mach number either side of 1. This 
is especially true of the flow in the immediate vicinity of the leading 
edge. Because of the nature of the flow field, the findings for the lower 
surface of the plate can be applied, within minor limitations, to the flow 
on the underside of any sharp-nosed airfoil with a flat lower surface . 

As in reference 3* the report is divided into two parts. Part I con- 
tains an outline of the general problem and its solution and a discussion 
of the final results. Part II supplies the mathematical details. A sum- 
mary of the notation is given in an appendix at the end of the report. 

The relaxation calculations for the lower surface of the plate were 
performed with unusual skill and diligence by Mrs. Marjorie Sill. The 
characteristics construction on the upper surface was programmed for the 
electronic computer by Meyer M. Resnikoff , who contributed many valuable 
ideas to this phase of the work. 

I - GENERAL PROBLEM AND FINAL RESULTS 


GENERAL PROBLEM 


Description of Flow Field 


The inviscid flow over an 
inclined plate at free-stream 
Mach number 1 has been described 
by Guderley in reference 1. A 
qualitative picture of the flow 
is given in sketch (a), which 
is essentially a reproduction 
of figure 1 of Guderley* s 
report. 

As shown in the sketch, 
there exists at angle of attack 
a large region of subsonic flow 
in front of and beneath the 
plate. This region is bounded 
by the lower surface DB of the 
plate, by sonic lines DO and BO 
springing from the leading and 
trailing edges, and by the free 
stream 0 at infinity. In this 
region, as in incompressible 
flow, the fluid forward of the 
plate is deflected upward. As 
a result, a stagnation point 
C occurs on the lower surface 
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of the plate near the leading edge. At this point the central (or "stag- 
nation") streamline branches to run both fore and aft along the plate. 
Sonic speed is reached on the surface of the plate precisely at the lead- 
ing and trailing edges D and B. The sonic line- from each of these edges 
leaves the plate at right angles to the lower surface. Downstream of each 
sonic line there occurs a supersonic expansion fan centered at the edge. 

Expansion fans of the type encountered here have been described in 
detail in references 3 and. 7 . In the limit of the edge itself, the fan 
approaches the classical Prandtl-Meyer flow. Away from the edge, however, 
the elementary expansion waves (or Mach waves), which are straight in the 
Rrandtl-Meyer case, now bend upstream toward the sonic line. Seme of 
these expansion waves, in fact, meet the sonic line, where they are 
reflected as compressions; others pass entirely downstream of the line. 

The wave that separates the two classes of expansion waves, the so-called 
"limiting wave" (see EO and AO in sketch (a)), approaches the sonic line 
at infinity* Obviously - and this is the essential point - the flow in 
the subsonic region must depend in part on conditions in the supersonic 
fan ahead of the limiting wave. It is, however, completely independent 
of conditions downstream of this wave. 

The details of the flow over the upper surface of the plate are not 
completely certain. Even in the absence of viscosity, a separated region 
must be expected adjacent to the leading edge. This follows from the 
fact that the angle of turn called for at the edge - l 80 ° in the present 
case - is greater than the 13 0° attainable by expansion to a vacuum. In 
any real fluid, of course, the pressure actually attained in the separated 
region will be fixed by viscous phenomena. Prom the standpoint of the 
present purely inviscid theory, however, this pressure is to be regarded 
as an assignable parameter (subject to one restriction to be mentioned 
later). Once the flow has separated from the plate, the central stream- 
line follows a path of constant pressure. To do this in the presence of 
the compression waves reflected from the sonic line, it must curve back 
again toward the plate. At some point the streamline will presumably 
reattach to the plate and be deflected to the direction of the upper sur- 
face by an oblique shock wave. Downstream of this reattachment wave, the 
flow will proceed at decreasing supersonic speed to the trailing edge, 
where a second shock wave will occur. The foregoing is, at least, a 
plausible description of the flow as it might be expected to exist on the 
upper surface. It is not inconceivable, however, that the shock system 
could be more complicated than that described here. 

The foregoing discussion implies one restriction that has not yet 
been stated. This is that the pressure assigned in the separated region 
must be not greater than that which exists where the limiting wave meets 
the leading edge. For pressures greater than this value, the l imi ting 
wave would have to start, not at the leading edgq, but at some point on 
the separated streamline. In this situation the part of the streamline 
between the leading edge and the foot of the limiting wave would be inter- 
dependent with the subsonic field. The problem would then be slightly 
different from that considered here, though the methods devised in this 
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report could aga 1 n be used to obtain a solution. Since the flow on the 
underside of the plate would he little affected in any event, this second 
possibility will be ignored, and the earlier situation will be presumed 
to exist. 

As should be apparent from this discussion, the analytical problem 
falls into two parts. To calculate the flow on the lower surface of the 
plate, a transonic boundary- value problem must be solved in a region con- 
sisting of the subsonic field and those portions of the two supersonic 
fans ahead of their limiting waves. Once conditions on the upper 3.1ml. ting 
wave are known, the calculation of the field above the plate is an initial- 
value problem in purely supersonic flow. In the first problem, the flow 
is everywhere irrotational . In the second problem, it is irrotational as 
far as the reattachment wave and rotational after that. 


Method of Analysis 

The method of analysis is described in detail in part II of the 
report. The mixed flow about the underside of the plate is found by 
numerical (i.e., finite-difference) solution of a boundary- value problem 
in the holograph plane. In this formulation of the problem, the free 
stream appears as a singularity - that is, an infinity - at a point on 
the sonic circle. This singularity consists of two superposed parts, one 
antisymmetric and similar to a doublet, the other symmetric and similar 
to a vortex. To calculate the flow, two separate boundary- value problems 
are solved, each containing one of the elementary singularities . The 
final result is then obtained by combining the two solutions to meet the 
condition that the central streamline branches at the stagnation point. 

As in all finite-difference work, the singularities must be incor- 
porated into the solutions analytically. This is complicated in the pres- 
ent case by the fact that simple closed expressions for the singularities 
are not available for the complete differential equation. A suitable 
method can be devised, however, based on the use of known results for the 
Tricomi equation. The method takes full account of the primary effect of 
the singularity. A small approximation is involved, however, in that 
singularities in certain of the higher- order derivatives (i.e., second 
order and above) are ignored in the numerical calculations. Experience 
indicates that the error from this source is probably smaller than that 
involved in the finite-difference process itself. It is sure to be 
negligible in the vicinity of the plate. 

The finite-difference solution of the boundary-value problems entails 
several elements of novelty. These are made necessary by the presence of 
the singularity and by the mixed nature of the flow field. In the sub- 
sonic region the finite-difference equations are established in what is 
essentially a polar coordinate system. Except for the inclusion of the 
singularity, the procedures here are more or less standard. In the super- 
sonic region the equations are set up in characteristic coordinates. The 
procedures here and at the transition from one coordinate system to the 
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other at th e sonic line are somewhat unusual. Solution of the system of 
finite-difference equations is carried out on desk calculators hy a com- 
bination of relaxation techniques and step-by-step procedures. Once the 
solution is obtained in the hodograph plane, transformation to the physical 
plane is a simple matter. 

The purely supersonic flow over the upper surface of the plate is 
found in the physical plane by a standard form of the method of character- 
istics. In the present work a completely numerical process was used, in 
contrast to the more usual semi graphical procedures . This was done for 
two reasons: (l) Tha high Mach numbers and correspondingly s m all Mach 

angles in the vicinity of the separated region made any semi graphical 
procedure of doubtful accuracy; and ( 2 ) the large number of points required 
in the characteristics net made the use of automatic computing machines 
mandatory. The numerical work was carried out on an electronic digital 
computer. As mentioned in the introduction, the calculations for this 
part of the problem have been carried only as far as the reattachment 
point on the upper surface. For the time being they have also been con- 
fined to only one value of the separation pressure. It is hoped that the 
work can be carried on to the trailing edge and repeated for other values 
of the separation pressure at a later date. 


RESUITS AMD DISCUSSION 


Calculations have been made, on the basis of the methods outlined, 
for an angle of attack of 13°. The work was carried out for a ratio of 
specific heats 7 of 1.4-05 instead of the more usual value of 7 / 5 . This 
was done because it was originally intended that the characteristics con- 
struction would be made by a semigraphical version of Guderley’s method, 
for which extensive tables on the basis of the former value are available 
(ref. 8 ). When this method was abandoned in favor of automatic numerical 
computation (which does not utilize the tables), the calculations for the 
lower surface with 7 = 1.4-05 were too far along to alter. The final 
results in the physical plane are shown in figures 1 and 2 and are dis- 
cussed in the following paragraphs. (For results in the hodograph plane, 
see part H of the report.) 


Flow Field 


Lower surface and leading edge .- The flow field about the lower sur- 
face and leading edge is shown in figures 1(a) through 1(d). To understand 
these figures, it will perhaps be best to follow them through consecutively 
in seme detail. 

The first figure (1(a)) shows the complete plate from the leading 
edge to the trailing edge. In this and the subsequent parts of figure 1, 
size is indicated by a scale of s/c, where s is any length in the field 
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and c is the chord of the plate. Streamlines are identified hy a value 
of [6(y/c) ] co) which denotes the vertical distance at infinity between the 
streamline in question and the central (or stagnation) streamline. Con- 
tours of constant fluid properties are identified by the value of the local 
Mach number M. 

The most striking thing about figure 1(a) is the extraordinarily 
sm all size of the region into which the large changes of flow around the 
leading edge are concentrated. This is typified by the distance from the 
stagnation point to the leading edge, along which the flow must accelerate 
from M = 0 to M = 1. Even at the fairly large angle of attack of 13°, 
this distance is only 0.0016 of the chord, which is too small to appear 
in the present figure. In incompressible flow about a lifting flat plate 
at the same angle of attack, the corresponding result as given by con- 
formal transformation (see, e.g., ref. 9) is 0.05 of the chord. Obviously, 
the stagnation point must move forward markedly as the free-stream Mach 
number increases in the subsonic range. When Mach number 1 is reached, 
it is, to the scale of the plate as a whole, practically coincident with 
the edge. The representation of the leading-edge flow by a singularity 
in which the stagnation point lies precisely at the edge is a well-known 
approximation in thin- airfoil theory. The present results suggest that 
the approximation should be even closer to the truth at transonic speeds 
(as, for example, in the work of Guderley (ref. 1)) than it is in the 
classical incompressible case. 

As with the displacement of the stagnation point, the details of the 
sonic line and limiting wave are not visible in figure 1(a). As nearly 
as can be seen from this figure, the sonic line meets the leading edge 
more or less directly from above, while the limiting wave comes into the 
edge slightly from the rear. This appearance is, of course, at variance 
with the qualitative description given earlier in connection with 
sketch (a) . If the results of figure 1(a) were all that were available, 
however, one might be inclined to accept the statements of the present 
paragraph as correct . 

That the true state of affairs is quite different begins to appear 
in figure 1(b) . This figure shows the flow over the forward portion of 
the plate to a scale 20 times that of figure 1(a) . The present figure 
also includes the separated streamline, which will be discussed later. 

Enlargement to the scale of figure 1(b) is sufficient to show the 
displacement of the stagnation point. A pronounced curvature of the sonic 
line is also visible here, and this line seems now to meet the leading edge 
from directly forward of the plate. The latter situation is, however, 
still different from that described in connection with sketch (a) . The 
limiting wave in figure 1(b) appears to meet the edge from the rear, but 
with greater slope and more curvature than were visible before. 
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To examine the situation at the leading edge in still greater detail, 
the results are replotted once more in figure 1 (c) , this time to a scale 
200 t ime s that of figure 1(a) . Here at last it appears that the sonic 
line does indeed approach the leading edge from the underside of the plate. 
The manne r in which the approach takes place, however, is still not clear. 

To see the latter details, one is forced ultimately to a plot such 
as figure 1(d) . This plot has a magnification 5000 times that of the 
original figure 1(a) . Because of this very large magnification, the stag- 
nation point in the present plot would lie approximately 50 inches off the 
page to the right. The trailing edge (whose position nevertheless deter- 
mines the characteristic length in the problem) would he slightly more 
than 1/2 mile away. Despite the very large scale of figure 1(d), the 
results that appear here can he specified with good accuracy. This is 
because the hodograph transformation has the property of greatly enlarging 
the region near the leading edge relative to the rest of the field. Thus, 
a perfectly reasonable mesh interval for the finite-difference scheme in 
the hodograph plane (see part H) can provide sufficient data to define 
the lines of figure 1(d) without difficulty. To obtain comparable accuracy 
in calculations in the physical plane an impractically small interval would 
be required in the vicinity of the leading edge. 

The contrast between the leading- edge flow as it appears in fig- 
ure 1(d) and as it appeared originally in figure 1(a) is obvious. In 
figure 1(d) it can be seen that the sonic line does in fact meet the lead- 
ing edge at right angles to the lower surface of the plate. (It can be 
shown that the curvature of the sonic line where it meets the plate is 
infinite.) The limiting wave in figure 1(d) comes into the edge from a 
direction slightly forward of the vertical. It thus approaches the plate 
much as the sonic line seemed to in the small plot of figure 1(a) . To 
the present scale it is also possible to see something of the way in which 
the influence of the leading edge is propagated in the supersonic region. 
This is shown by the Mach lines included in the figure. The reader may 
find it interesting to compare these results with a similar plot of the 
Mach lines and streamlines in the classical Prandtl- Meyer flow (see, 
especially,^. 278 of ref. 10 and p. 171 of ref. 11). 

Separated streamline .- The separated streamline is shown in fig- 
ures 1(b) through 1(d) for the single case that has been calculated. In 
this case the pressure p g on- the separated streamline is given by 
Ps/p-too = O.OOO738, where p^ is the total pressure in the free stream. 

This is precisely the value of the pressure ratio that exists where the 
limi ting wave meets the leading edge - that is, we have assumed that the 
flow remains attached to the leading edge until the limiting wave is 
reached and then immediately separates (see "Description of Flow Field"). 
The corresponding Mach number M s on the separated streamline is 5*88, 
and the angle at which the streamline leaves the plate is 83-1/2° relative 
to the free-stream direction. 
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As can "be seen from figure 1(b), the over- all dimensions of the 
separated region are of the same order of magnitude as the distance from 
the leading edge to the stagnation point. The length of limiting wave that 
is needed to calculate the separated streamline is marked off in fig- 
ure 1(h) by the point P. This is the point from which the downgoing 
Mach wave to the reattachment point leaves the limiting wave. The angle 
through which the flow must be deflected at the reattachment point turns 
out to be 19°. This is well below the maximum deflection that can be 
attained through an oblique shock wave in an inviscid flow at the assumed 
Mach number of 5*88. Except for these observations, little can be said 
about the separated flow until results become available for other values 
of the pressure in the separated region. 


Pressure Distribution 


The distribution of pressure coefficient Cp on the lower surface 
of the plate is shown in figure 2. Results for the complete surface are 
given in figure 2(a) and for the region between the leading edge and stag- 
nation point in figure 2(b). Included for comparison in figure 2(a) are 
three additional sets of results obtained an follows: 

(a) Directly from Guderley’s smal 1 -disturbance analysis for M ra = 1 
(ref. 1), 

(b) From Guderley’s small-disturbance analysis with the approximate 
relationship between his holograph variable r\ and the dimen- 
sionless speed w replaced by an exact relationship (see 

eq. (10a) of part II; this procedure wan suggested and used by 
Guderley for a different problem in ref. 12), 

(c) From the classical methods of conformal transformation for 
incompressible flow (see, e.g., ref. 9) • 

It can be seen from figure 2(a) that Guderley’s work of reference 1 
somewhat overestimates the pressure over the entire chord. The average 
value of Cp given by his analysis is 0.59; that of the present work is 
O. 525 . According to Guderley’s calculations, the center of pressure is 
at 35 percent of the chord; the present work puts it at 37-9* Considering 
the fairly large angle of attack, however, the quantitative accuracy of 
the small-disturbance analysis is remarkably good. The only major error 
in the small-disturbance analysis occurs near the leading edge where the 
approximate theory represents the stagnation point by an infinity at the 
edge. The present results, by contrast, rise to the stagnation Cp of 
1.275 very close to the edge and then drop to zero at the edge itself 
(see fig. 2(b)). The effect of this qualitative difference on the inte- 
grated force on the surface is obviously small. 

According to figure 2(a) , application of a correction to the small- 
disturbance results as outlined in (b) above eliminates most of the 
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nume rical error relative to the present findings. The discrepancies that 
do remain are, for the most part, within the limits of accuracy to which 
one can read numerical results from the graphs of reference 1. As regards 
the leading-edge flow, the corrected curve gives a true stagnation value 
of 1.275 "but places the stagnation point still on the edge. For almost 
all purposes, therefore, there appears to "be no need to go "beyond the 
corrected small-disturbance analysis. Only if the details of the flow in 
the immediate vicinity of the leading edge are of importance is it neces- 
sary to resort to more accurate work. 

The results for incompressible flow have "been included in figure 2(a) 
for general interest. They illustrate the relatively great displacement 
of the stagnation point (Cp = 1.0) that occurs in the incompressible case. 


CONCLUDING REMARKS 


To what extent the results found here will apply to an actual airfoil 
in a real fluid is an open question. The only experimental data bearing 
on the problem are those of Wood (ref. 13) , who made interferometric 
studies of the flow at high subsonic speeds over a thin wedge at an angle 
of attack of the lower surface of 11. 5°* Wood's results for the lines of 
constant Mach number (see, in particular, his figure 3(e) for Moo = 0.89^) 
show a distinct resemblance to the flow field of figure 1(a) . Wood remarks 
on the basis of his results that the stagnation point is "very close to 
the leading edge," though precisely how close his measurements do not show. 
At moderate distances above the leading edge, the sonic line found by Wood 
looks reasonably like that of figure 1(a) if regard is had for differences 
in free-stream Mach number. As to conditions in the immediate vicinity 
of the edge. Wood takes his results to "indicate strongly that the sonic 
line starts from the upper surface," though "refraction of the light in 
the large density gradient . . . precludes following the contours all the 
way to the surface." Because of the latter circumstance, the interfer- 
ometer, in effect, views the field from essentially the scale of fig- 
ure l(a) (or, at best, something a slight bit larger) . In view of the 
minute size of the leading- edge flow as calculated in the present work. 
Wood's inference must therefore be taken as tentative. As factors that 
might move the sonic line from the lower surface. Wood cites the viscosity 
of the fluid and the nonzero thickness of the real leading edge. Viscos- 
ity will undoubtedly be of some influence, though the large negative 
gradient of pressure near the leading edge would tend to minimize this 
effect. (The pressure in the separated region will, of course, be fixed 
by viscous phenomena, but these need not influence conditions on the lower 
surface . ) The thickness of the leading edge is likely to be of more 
importance, since any real edge, though "sharp" in the ordinary sense, 
will have an enormous thickness when viewed to the scale of figure 1(d) . 

It may well prove, however, that any changes from the the oreti ceil pattern, 
even from this latter source, are of only local effect and that the rela- 
tively distant field will be correctly given by the theory. Certainly 
there are many questions to be answered, and the remarks of this paragraph 
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are to "be considered merely as miscellaneous thoughts on the subject. A 
challenging field, exists here for experimental research. 


II - DETAILS OF ANALYSIS 


An outline of the method of analysis has been given in part I of the 
report. The material that follows is concerned with the mathematical 
details. In the initial section the analytical basis is given for the 
calculation of the mixed flow over the lower surface in the hodograph 
plane. The main interest here is in the treatment of the free- stream 
singularity. The second section is concerned with the numerical solution 
of the hodograph problem. The procedures developed here, particularly as 
regards the finite-difference equations in the supersonic region near 
the sonic line, may have application in other problems of mixed flow. 

The third and final section deals briefly with the characteristics con- 
struction used to find the purely supersonic flow over the separated region 
on the upper surface. 


ANALYTICAL BASIS FOR CALCULATION OF FLOW OVER LOWER SURFACE 
Boundary- Value Problem in Hodograph Plane 


. A representation of the flow over the lower surface of the plate in 
the physical plane has been given in part I (see sketch (a)). The corre- 
sponding picture in the hodograph plane is shown in sketch (b) . The 
variables here are the dimension- 
less speed w (i.e., the ordinary 
speed V made dimensionless 
through division by the critical 
speed a*) and the flow inclination 
0 (measured relative to the free- 
stream direction) . Corresponding 
points in sketches (a) and (b) are 
noted by the same symbols. 

In the hodograph plane the 
sonic speed appears as the circle 
w = 1. The free stream appears 
on this circle as the point 0 
located at 0=0. All of the 
streamlines issue from this 
point, which is therefore sin- 
gular. The lower surface of 
the plate is represented by a 
straight line passing through 
the stagnation point C and 
inclined at an angle of attack 



Sketch (b) 
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a to the free- stream, direction. This line meets the sonic circle at 
points D and B, which correspond to the upstream side of the leading and 
trailing edges, respectively. The leading edge itself appears as the 
clockwise characteristic (i.e., epicycloid) through Dj the trailing edge 
as the counterclockwise characteristic through B. The limiting Mach 
waves appear in the holograph as the. two characteristics passing through 
the free -stream point 0. The point (E or A) at which each of these char- 
acteristics intersects the image of the corresponding edge fixes the 
extent of the edge that is of importance in determining the flow over the 
lower surface. The stagnation streamline rims in the holograph from the 
free- stream point 0 to the stagnation point C, where it "branches. The 
streamlines that lie above the stagnation streamline in the holograph 
plane pass above the plate in the physical plane. Those that lie below 
this streamline pass below the plate. 

To calculate the flow over the lower surface of the plate, a boundary- 
value problem must be solved in the region QABCBEO. If the stream func- 
tion Y(w,0) is taken as the dependent variable, the differential equation 
to be satisfied is the usual linear equation (see, e.g., ref. lk, p. lk7) 


W2f w + W 


y-3 

7 + 1 


w 2 


= 0 


( 1 ) 


where differentiation is indicated by the subscript notation. This equa- 
tion is elliptic for w < 1 and hyperbolic for w > 1. The characteristics 
that exist for w > 1 are given by 


or 



(2a) 



where ©w=i is the value of 0 
sonic circle. 


at ■which the characteristic meets the 
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The 'boundary conditions to "be satisfied in the problem are as fol- 
lows: 

1. The value of I is constant - say 0 - on the boundary ABODE. 

2. The function Y has the proper singular behavior at the free- 
stream point 0. 

3. The streamline I = 0 branches at the stagnation point C. 

In order to represent the free stream, the singularity at 0 must have 
the following properties (cf. ref. 15 ): 

(a) The function Y takes on all values from +00 to -00 at the 
point 0. 

(b) Wo singularities propagate along the limiting characteristics 
EO and AO. (This follows from the fact that no singularities 
exist in the boundary conditions at E and A.) 

(c) The limiting characteristics transform into the finite part of 
the physical plane. 

(d) The flow maps onto a single sheet in the physical plane. 

Simple and useful expressions satisfying these requirements for the small- 
disturbance equivalent of equation (l) (i.e. , the Tricomi equation) have 
been given by Guderley (refs. 1, 15, and 1 6 ) and Frank! (ref. 2). These 
expressions are the Siam of two parts: (l) a doublet-like singularity 

antisymmetric in 0 and (2) a vortex-like singularity symmetric in 0. 
Correspondingly simple expressions for the exact equation (l) have not 
been given, though Frank! (ref. 2) has obtained an infinite series of 
Chaplygin solutions that satisfies the given requirements. The present 
work will utilize the known singularities for the Tricomi aquation (which 
predominate in any event) plus certain correction terms sufficient to 
account for a significant part of the difference between the approximate 
and exact results (see following section). 

In view of boundary condition 1, boundary condition 3 can be met by 
requiring that 


¥ W (O,0) = 0 


(3) 


for any (and hence all) 0 between -a and jt - a. To satisfy this condi- 
tion in the present work we shall write ¥ in the form 


(M 


Y = f 8 - + c? 8 
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where S 8, and f 8 are functions containing, respectively, the antisymmetric 
and symmetric parts of the free-stream singularity and C is a constant. 
The functions S 8, and will "be required to 'satisfy individually the 
"boundary condition 1 "but not the "boundary condition 3* The combined func- 
tion ¥ will -then constitute a solution of the complete problem - that is, 
will satisfy all the boundary conditions including 3 - provided C is 
evaluated such that 


?w a (o,e) 

*w S (o ,0) 


(5) 


The task now is to solve the individual boundary-value problems for 
'S a and I s . In the end this will be done by numerical (i.e. , finite- 
difference) means . First, however, certain analytical preliminaries are 
required, particularly with reference to the free-stream singularity. 


Transformation of Differential Equation 


To utilize the known singularities for the Tricomi equation, it is 
necessary to put equation (l) into a form closer to that of Tricomi. 

This can "be done "by a transformation of the type used "by Guderley in the 
case of the equation for the Legendre potential (refs. 12 and 16): 


q = tj(w) 
1 


* = 


g(v) 


Y 


(6) 


Expressions for the functions q and g appearing here are found "by sub- 
stituting equations (6) into equation (l) and requiring that the resulting 
partial differential equation for ijr agree in its derivatives with the 
Tricomi equation - specifically, that it have the form 


*qT] “ W e9 + = 0 


(7) 


This procedure leads to the following ordinary differential equations for 
q and g: 


w 2 - 1 


a - ^ v=, 


7+1 




(8) 
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(9) 


where the primes denote differentiation of a function with respect to its 
argument. By integration of equation (8 ) , tj is found as 

/ "> 2/3 



(10a) 

for w < 1 and 



(10b) 


for w > 1. Here T] has been taken equal to zero at w = 1 to satisfy 
the requirement that equation (7) change type at this point. (The simi- 
larity between equations (10b) and (2b) should be noted.) Pram equa- 
tion (9), g is found as 


7+2 


. 1/4 


g(w) = 2 


1/6 


6(7-1) 

(¥) fe) (* - 


7-1 
7+ 1 


w 2 ^ 


7 + 1 
4(7-1) 


( 11 ) 


where the constant of integration has been chosen (arbitrarily) to make 
g = 1 at w = 1. Plots of T] and g as functions of w for 7 = IA05 
are shown in figures 3 and 4, and tabular values are given in table I. 
The function t] is positive at supersonic speeds (w > 1) and negative 
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at subsonic speeds (w < 1) . It goes to negative infinity as w — > 0 and 

has a finite value (r) fflOT ) at the maximum speed Vma-v - *J(7 + 1)/ (7 “ 1) • 

The function g goes to positive infinity as w — > 0 and to zero at 
"Wmax- 

The function b(t]) in equation (7) is given Implicitly "by. the equation 






£. 

s 



( 12 ) 


where g’/g g n /g are to "be evaluated with the add of equations ( 9 ) 
and (8) . A graph of b as a function of i) is given in figure 5; 
tabular values are included in table I. The value of b is always nega- 
tive; it approaches negative infinity as t\ — > ^ and zero as tj — > -<». 

For the subsequent work it is also necessary to know b as a power series 
in tj about tj = 0. The required expansion is 


ToC-q) = b Q + biti + b 2 tl 2 + . . • 


(13) 


where 


b Q = -0.9181 
b x = -1.3164 


b 2 = -1.0741 


e 



When transformed 
into the tj,0 plane, 
the holograph boundaries 
appear as shown in 
sketch (c). The char- 
acteristics of the dif- 
ferential equation (7 ) , 
which define the bound- 
aries for T) > 0, are 
given by (see eqs. (2b) 
and (10b)) 


Sketch (c) 



MCA TN 3723 


17 


0 = 0rj=o ± f ^ 3/2 


(HO 


These are the same, of course, as for the Tricomi equation. 

To carry out the solution for tj > 0 we shall also require equa- 
tion (7) in characteristic coordinates. These coordinates will he denoted 
here by r and l (corresponding to the right- and left-hand characteristics 
relative to the direction of positive T|). They are defined by the 
equations 


or equivalently 



(15a) 

(15b) 

(16a) 

(16b) 


When transformed into these coordinates, equation (7) becomes 


[6(r- l)f /3 -t rl + z) (* r - * x )_ 


+ bi|r = 0 


(17) 


Free-Stream Singularity 


Method of treatment .- The finite-difference solution of problems 
involving a strong singularity (i.e., an infinity or discontinuity in the 
unknown or an infinity in its first derivatives) has been discussed by 
various writers (see, e.g., ref. 17). The usual procedure in linear prob- 
lems is to write the dependent variable as the sum of two parts, each of 
which is individually a solution of the differential equation. The first 
part is a term of known analytical form containing the required singular- 
ity. The second part, free of singularity, is to be determined in accord- 
ance with the remaining, nonsingular boundary conditions. This part, being 
regular, can be found by numerical means. This procedure, unfortunately. 
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cannot "be followed in tbs present case , since suitable analytical expres- 
sions for the singularity for the complete equation (7) are not known. 
Instead the procedure here will be to utilize the known singularities for 
the Tricomi equation plus certain analytical correction terms. These 
terms will reduce the order of the singularity in the remaining unknown 
to the point where it may safely be ignored in the numerical process. 

The sonic singularities for the Tricomi equation - that is, for the 
equation 


T(t) = " Tplr 00 = 0 (18) 

have been discussed at length by Guderley in reference 16. Following 
Guderley’s notation, we shall express these solutions in the form 

+ = jn^fjC^jn) 


where £ is a new variable defined by 

r - " 2 

V 3 

The function f j is a solution of a hypergeometric equation whose coef- 
ficients involve n as a parameter. Expressions for fj are listed in 
reference 16 in terms of hypergeometric functions. As shown in refer- 
ence 16 , the values of n required to represent a uniform free stream 
are -5/2 for the antisymmetric singularity and -1 for the symmetric 
singularity. 

To utilize the solutions (19) in the present case, we begin by 
rewriting equation (7) as 


(19) 


( 20 ) 


T(tJt) = " (130 + "bin + ^a 1 ! 2 + . . .)t 


( 21 ) 


For reasons that will appear directly, the solution of the boundary- value 
problem (for either ty 8, or' i]r s ) is then taken in the form 


* 



+ X 


( 22 ) 
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where A is any convenient constant. The ^ that appear here are 

singular correction functions 2 of the form ^ = iTjf^g^C^) and X is 

added to satisfy the nonsingular houndary conditions. In contrast to the 
us u al procedure outlined above, none of the terms in equation (22) are 
individually solutions of the differential equation (21) . To make? clear 
the reasons^? or the foregoing choice - anfl at the same t ime derive equa- 
tions for i | r and X - we substitute equation (22) into equation (21) and 
arrange the result according to the following scheme: 


- [T(X) + bX] 


f 

+ T(y c ) 


* 0 * 

+ t(?j) 


+ 

+ T(? 2 ) 

H 

+ + i3o9 0 

+ t(? 3 ) 


+ ksn 3 ? + 

+ « • . 



+ T(+ m ) 


+ VP? + VvP'^o + VsHi + . . . + b Q v 2 

+ + Wf 2 ?! + . . . + b lT) £ m _ 2 + b 0 $ arl 

+ + VPto + ViHj + . . . + WV;, + + k 0 v 



(23) 


The idea here is that we determine $ 0 , ij^, etc., such that T($ 0 ) = -b 0 ijr, 
T($ x ) = -bjTpjf, etc. Since, in general, T^) = |tj | a±- ^Ix t (g) if 

= h | ni gi(£), ^ since we have ijr = |n pfj, it follows that the correc- 
tion functions must proceed like = [tj | n+2 g Q , $ x = |rj| Iws g ;L , etc. The 

scheme of equation (21) follows accordingly. (Note in particular that for 
1 > 2 the previous correction functions must be included together with ijr 
on the right-hand side of the equation. ) It thus appears that the expres- 
sion (22) will constitute a solution of the differential equation (21), 
provided the general correction $1 and the function X satisfy the fol- 
lowing nonhomogeneous differential equations: 


Correction functions of this type have also been used by Frankl in 
reference 2. 
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for 

T(X) 


T(*i) = - (b^ 1 * + b.£_ 2 T]i- 2 £o 

+ ' b i..3 T l i " 3 $i 

i = 0, 1, . . . , m and 


+ bX = - ACO^ai pH. + b^Tfi +2 

+ - . • )$ 

+ (bm- lif 1 " 1 + bjnif 1 

+ . - - )to 

+ (bm-2Ti m_2 + bm-iTf 1 ” 1 

+...)?! 

+ 


+ (b 0 + bin 

+ . . . ) . 

+ (bo + "ball 

+ . . • )t m ] 


+ "boti-a) 

(2k) 


( 25 ) 


It is apparent that the correction functions like the original func- 

tion ijr, are singular at the origin. The singularity, however, becomes 
weaker (i.e., makes its first appearance in progressively higher-order 
derivatives) as i increases. The function X is also singular at the 
origin, the singularity here being one order weaker than that in ijr m . 

Equations (2k) and (25) form the basis for the present work. As 
will be seen, equation (2k) can be solved analytically for as large a 
value of i as will be required. Solution for arbitrary i would be 
very difficult. (This is the reason for not taking m = co and thus 
reducing X to a regular function.) Equation ( 25 ) will be solved numeri- 
cally, subject to appropriate boundary conditions, by ignoring the singu- 
larity in X at the origin and applying the finite-difference techniques 
ordinarily used for a regular function. If the singularity is sufficiently 
weak, the resulting error in the over-all solution should be small. To 
make a numerical solution possible at all, a mi nimum requirement is that 
the right-hand side of equation ( 25 ) must be finite and single-valued at 
the origin. Since ijr = |n ^fj, it can be seen that for this requirement 
to be satisfied we must choose m such that m + l + n > 0 or that 
m > -(n + l). For the antisymmetric singularity (n = - 5 / 2 ) we shall 
therefore take m = 2. For the symmetric singularity (n = -l) , m will 
be taken equal to 1. In the antisymmetric case the singularity that then 
remains in X first appears in the derivative Xqq (but not in or 

X^q) . In the symmetric case it first appears in the third- order deriva- 
tives. Ignoring singularities of this order in the solution for X may 
be expected to cause negligible error, especially in the vicinity of the 
plate . 


Solution for correction functions .- According to the foregoing, 

' expressions for the correction functions ^ are needed for i = 0, 1, 2. 
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These are found "by solving equation (2^) , -which, can also he written 

i -2 

$i - ^i 00 = “ ^i 1 ! 1 * + Yi ^-a-t 1 ! 1 "*" 1 ^ (26) 

t=o 


To carry out the solution we introduce a function tjr that satisfies 
the Tricomi equation 

*in - ^00 = 0 ( 27) 


and has the property that 

lee = 1 (28) 

(Expressions for ijr will he given later.) 3 From equations (27) and (28) 
we also have the relations 


♦ip, = *1*00 - ^ ( 2 9 a ) 

and 

\ee = *00^ “ *n ( 2 ^ 


The solution of equation (26) will now he assumed in the form 

$i = PiT] 1-3 ^ + (30) 


where p^ and q^ are constants whose values are to he found. By differ- 
entiating this expression and making use of relations ( 29 ) , one obtains 


s This approach was suggested hy Gottfried Guderley, to whom the 
authors are much indebted. 
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= P ± [(i- l)(i“ + 2(i - ih 1 " 2 * + ] + 


’•m 


q. ± [i(l- Dn 1 - 2 ^ + (21 + l)^ + Ti i+1 tT]] 


( 31 a) 


$iee - Pi T i i “ ll J r e0 + 


(31b) 


Substitution of expressions (31) and (30) into equation (26) then gives, 
in view of equation (27), 


q i (2i + l)i] i iir + P ± (i- l)(i - 2)Tj i-3 t + (2p± + iq ± )(i - 


i-2 


i-2 


= -hiql$ - ti 1 -- 3 * Y, \-2-iPt ~ X Va-A 

t=o t=o 


For this equation to he satisfied for all r\ we must have 


q i (2i + 1) = -hi 


Pi (i - l)(i - 2) = -X Va-t^t > 


(32) 


t=o 


i-2 


(2p ± + iqi)(i - 1) = - X b i- 


-t^ 


t=o 


These are three equations for the two unkn owns Pi and qi, so that in 
general a solution cannot he found. For i = 0 , 1 , 2 , however, certain 
of the equations either coincide or disappear, and a solution can he 
obtained as follows: 
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For i = 0 equations (32) reduce to 

<3o = -l3 ° 

"2P 0 = 0 

-2p 0 = 0 

which have the solution p Q = 0, ^ = -b 0 . 

For i = 1 equations (32) become 

3 ( 1 ! = 

p x . 0 = 0 

(2P! + qj .0 = 0 

We thus have p x arbitrary - say 0 - and q 1 = -bi/3. 4 
For i = 2 equations (32) become 

5<k> ■ _ ' b 2 

p 2 . 0 = -boP 0 

2p 2 + 2qg = -boqQ 

Since p Q = 0 and q^, = -b 0 , these have the solution q a = -b 2 /5 and 

P 2 = (^o 2 /2) + (ba/5). 

By using these results in equation (30), we thus obtain the required 
correction functions as 



4 It is to be expected that p x would be arbitrary, since for i = 1 
the first term of expression (27) reduces to p x ^r, and this term will, by 
virtue of equation (2^) , disappear when substituted into the left-hand 
side of equation (23) , irrespective of the value of p 1 . 
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If it were required, a solution for \j? 3 ^could be found by taking a suit- 
able value of p x / 0 . A solution for $ 4 (and higher) in the form of 
equation (27) is apparently impossible. 

To apply equations (33) it remains to write expressions for 
and i}^ . With the aid of the equations of reference 1 6 (pp. 58 and 3^-35) 
it can be shown that, if f = |q | fj(£;n) as given by equation (19) , then 
an expression for ijr that satisfies equations (27) and (28) is 


_ , ,n+3 

\Jr = K|ti | fi(^n + 3) 


( 3k ) 


where 


(4n + l)(4n + 7) 

(n + 2)(n + 3)(2n + l)(2n + 3) 


By differentiating equation (34) and having recourse again to the equations 
of reference l6, it can also be shown that 


^ = K|ri| n4S (sign q)^(n+3)fj({;;n+3) - |£| (sign 0)fi(£jn + |^J 

(35) 

where (sign q) or (sign 0) are quantities with absolute value 1 and the 
sign of q or 0. 

Final relations .- By using equations (19), (33), (3*0, and (35) in 
conjunction with equations (22) and (25) , we can now write the final 
relations to be used in solving the boundary- value problems for \|r a and i|r 8 . 

Problem for i|r a : In this case the solution must contain the anti- 

symmetric singularity, in which case n = -5/2 and m = 2 . Equation ( 22 ) , 
written for thus becomes, after substitution from equations (33), 

¥*■ = A a [V+ (^- + qt a - (bo + y q + y 4 2 ) \ ] + X a 


(36) 
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The differential equation for the unknown X a is found "by substituting 
equations (33) into equation (25) and can be written 

Xtjti “ nXgg + tx£L = “ A a -|[b - (b 0 +b 1 T] + b 2 T] 2 )3t a + b^|- + t)f a - 

K*o + “ 1t 1 + y I 2 )' b o 2 ]?* } (37) 

The formulas to be used in computing if? , if , and ipf 1 , as obtained from 
equations (19), (3 1 *-), and (35), are 



^ = 'X H 1/2 ( s ign n) 0 " 3 l^| 1/2 ( sign _1 ) 


These quantities are al 1 antisymmetric with respect to 9 , so that 
computations need be carried through only for positive 0. 

Problem for ijr 8 : Here the solution contains the symmetric singular- 

ity, in which case n = -1 and m = 1. The equation for ijr s is therefore 

* 8 = A s |y - (b 0 + ^ ti) f ® ] + X s (39) 


and the differential equation for X s is 

x ^rj “ ^*00 + hY - a = " AS ^ b “ (ho+hifl)]^ 8 - b ^b 0 + y T]) f 8 j- 


m 
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The formulas 


~S — Q 

for ^ an d ijr^ 


are 


-1) 

= t i[9T i (6s 2) - 6|^j 1/2 (sign 0)f^J |)] 


(41) 


These q uant ities are "both symmetric with respect to 0. 

Formulas for the fj(£jn) that appear in equations ( 38 ) and (hi) are 
given by Guderley in reference l6 in terms of hypergeometric series. For 
the values of n required here, however, the results can be put in closed 
form analogous to that given by Frank! for n = - 5/2 in reference 2. If 
a new variable s is introduced according to the definition 


s 




(42) 


the pertinent formulas are 


<<• - i) 

*1 (is -i)- 

f ^ ! s ) 


*i(£;2) 


-I! k - 1 i 6 [(l + B ) x ■ s)S • (5 ' B ) (1 + 8)3 ] 

-| |«« - l|®[d - a) 3 + (1 + a) 3 ] 

- - \ | a 2 - 1| °[(1 - .)* - (1 + .)*] 

- - f |a a - lf 5 [(| + «)l - => 3 + (5 - + =) 4 ] 






J 


(43) 
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Equations for Transformation to Physical Plane 


After the boundary-value problems for tJt 8, and. t s are solved., W is 
found, from equations (6) and. (4) , and. the solution is transformed to the 
physical plane. If x and y form a right-hand coordinate system with 
x in the free- stream direction, the transformation equations can be put 
in the form (cf. ref. l4, pp. 146-147) 


where 


dx = Xtfdw + XqcL© 
dy = y v dw + y 0 d0 


x w — 


r o m 1-w 2 vj 

sin 0f v + cos 0 — 

L 1 y ~ 1 J 


1 - - — - w 2 
7 + 1 


x 0 - 


pwa* 


cos 0? w - sin 0 T 0 


y w pwa» 


r 1 - w 2 Ygi 

cos 9 I w sin 0 

L w , 7- 1 - W J 


1 - 1 — - w 2 
7 + 1 


r + etj 




In the present work al 1 distances will be made dimensionless in terms 
of the airfoil chord, which can be written 


X B ~ X D 
cos a 
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where xp and X 3 are the x coordinates of the leading and trailing 
edges (measured. relative to any convenient origin). By. applying equa- 
tions (44) through (46) along the line DCB, for which d0 = 0 (see 
sketch (h) ) , one can write 


X B " X D = 



m 


where 


I s - 



1 - V s *0 



y- 1 
y + i 



w 2 


d(±w) 


(49) 


It is understood from the limits that the integration 1 b taken from D 
to B along the image of the plate. The minus sign in the differential 
d(±w) is to he used from D to C and the plus Bign from C to B. Combi- 
nation of equations (47) and (48) then gives 


c 



(50) 


from which we can write 



j (a^x^w + a*x 9 de) 

> 

J (a^y w dv + a^y 0 d0) 


(51) 


Equations (51)# (49)# and (45) are the basis for the numerical transfor- 
mation from the holograph plane to the physical plane. 

To transform the characteristic lines in the supersonic region, it 
is convenient to rewrite the transformation equations in the coordinates 
r and l. This can be done on the basis of equations (15) and (8). The 
results are 
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where 


and 


d(f) = j (a*dr + a^dz) 
*(f) = J (a*y r dr + a*y z dz) 


( 52 ) 


= 


Pt 

pwa* 


(3 cos 0 - sin 0)¥f 




X Z - - 


(3 cos 0 + sin 0)^ 


Pt 

pwa* 

y„ = (P sin 9 + cos 0)¥ 

r pwa* ' r 




p t 

y Z pwa* 


(3 sin 0 - cos 0 )y 7 


J 


(53) 


1/2 


P = 


W 2 - 1 


a-^w 2 


7 + 1 


NUMERICAL SOLUTION FOR FLOW OVER LOWER SURFACE 
General Scheme 


The numerical solution of the boundary- value problems for ijr 8, and i]r s 
will be based on equations ( 36 ), (37) an d (39) j (^0). To expl a i n the 
general scheme, it is convenient to introduce the linear differential 
operator L^(F) defined equivalently by 
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^iyri “ + "KF 


L S (F) s 


[6(r-Z)] 2/3 [-F rZ + -- -[ -- (P r -F z )] + hF 


(5^a) 

(54b) 


In this notation the differential equations ( 7 ) and ( 17 ) become 


Ifc(+) = 0 (55) 

It is also convenient for present purposes to write the solutions (36) 
and (39) in the form 


t = AM + X (56) 


where M is a shorthand notation for the quantity appearing in brackets 
in the earlier representations. By substituting this expression into 
equation (55) , we can then write the differential equations (37) and (4o) 
in the form 


I^(X) = -AL d (M) ( 57 ) 

where L^(M) is identical to the quantity in braces on the right-hand 
side of the earlier equations. 

To carry out a solution for the unknown X, we begin in the us ual 
fashion by replacing the differential expression on the left-hand side 
of equation (57) by a corresponding difference expression. This is done 
as usual by covering the region of solution with a suitable mesh and 
approxi ma ting the differential operator L^(F) by an appropriate difference 
operator at each mesh point. This difference operator, denoted by L a (F), 
will appear in the form 


L A (F ) - L a (F 0 , Fj_, . . . , F n ) 

where the Fj (j - 0 , 1 , . . . , n) denote the values of F at the point 
in question and at n adjacent points. In this manner the differential 
equation (57) is replaced at each mesh point by a difference equation of 
the type 
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L a (X) = -AIfc(M) ( 58 ) 

So far the procedure is conventional. At this point one could pro- 
ceed, also in conventional fashion, to solve the difference equation (58) 
for X. The required values of X on the boundary would follow directly 
from equation (56) and the known condition that ijr = 0 on the boundary. 
Following the determination of X at all mesh points, the corresponding 
values of ijr would then be calculated, again by means of equation (56) . 

Considerable simplification over the foregoing method can be obtained, 
however, by using a method due to Woods (ref. 17 ). According to this 
method, instead of solving for X, one uses equation (56) immediately to 
replace Xj for all mesh points in equation (58) by the equivalent 
\jrj - AMj. (The special case of equations involving the origin, at which 
Mj « 00, will be treated later.) In this manner one obtains a system of 
difference equations of the form 

L A (*) = -A[Ifc(M) - L A (M)) ( 59 ) 

This system can be solved directly for t|t. This method, though completely 
equivalent mathematically to that outlined in the preceding paragraph, has 
a great advantage for the present work. This stems from the fact that the 
right-hand side of equation (59) , which is proportional to the difference 
between Ljj and L A both operating on M, tends to zero as the distance 
from the origin increases. At some distance (depending upon the value of 
A, the size of the mesh interval, and the accuracy desired in the work) 
the right-hand side will, in fact, become negligible. Beyond this dis- 
tance equation (59) reduces for all practical purposes to 


L A (t) = 0 (60) 

(which is the same as would have been obtained if the finite-difference 
approximation had been introduced into equation ( 55 ) ) • This means that 
the values of L^(M) and M need be calculated at only a small percentage 
of the mesh points, instead of at all points as would be the case if the 
method based on equation (58) were used. Furthermore, the advantages to 
the computer of working with the familiar variable i{r, with boundary con- 
ditions directly in this variable, are considerable. 

The foregoing procedure must be modified when one of the X j in 
the difference equation (58) - say X^ - is located at the origin. This 
situation will prevail in the difference equations at the origin itself 
and at points adjacent thereto. In these equations X^ cannot be replaced 
by \Jfj_ - AMjj. since is infinite. It must therefore remain in the 

equation, although the other Xj can be replaced as before. Even in such 
cases, however, the general procedure can be formally retained if we look 
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upon X]j as a fictitious value of ijr - say ~ and take Mjj. as arbi- 
trarily zero. If this convention is used in the expressions to be given 
below for L^(F) , equation ( 59 ) will automatically provide the correct 
difference equations at the origin and adjacent points. The value of ^ 
that is computed, however, will have no real significance. 


Details of Finite-Difference Equations 


To apply equation ( 59 ) , it is necessary to obtain expressions for 
the difference operator L a (F) . This will be done differently in three 
different regions of the hodograph, as illustrated schematically in 
sketch (d) . (Here only positive 9 will be discussed. The procedures 



Sketch (d) 

for negative 9 will be obvious.) For purposes of the present discussion 
tw is used to denote some convenient negative multiple of A, where A 
is the interval of the finite-difference mesh in the t \,9 variables; 
is defined according to the relation + A. With this notation, 

the three regions in the hodograph and the procedures used in each are 
briefly as follows (details will be given later) : 

Tjj < ti Is 0: In this region the field is covered with a square mesh 

of basic interval A in the t] and 9 directions. Adjacent to the upper 
boundary the interval is adjusted (made greater than A) so that the 
terminal mesh points lie on the boundary. Throughout the region the dif- 
ference operator is found by approximation to the differential operator 
in the form given by equation ( 5 ^a) . 

q > 0: In this region the mesh is formed by lines of constant 
r and l as shown in the sketch ( see also eqs . ( 15 ) ) . It is convenient 
here to include a line of constant r beginning on the sonic line a dis- 
tance A above the uppermost line in the subsonic field. At all points 



NACA TN 3723 


33 


(except those next to the sonic line) the difference operator is obtained 
by approximating the differential operator in the form given by 
equation ( 5 ^b ) . 5 


tj < Tij-j-: Here a reversion is made from q to the original vari- 

able w. This is done to avoid difficulties arising from the fact that q 
extends to -00. A rectangular mesh is employed with basic interval A in 
the 9 direction and A* in the w direction. At the right-hand side 
of the region the intervals in w are adjusted so that vertical mesh 
lines fall at locations Wj and Vjj corresponding to rjj and tj-j-j-. At 
all points the finite-difference equations are based not on equation ( 59 ) 
but on the differential equation (l) for 1 . For this to be permissible, 
rjj must be taken far enough to the left that equation ( 59 ) > if used, 
would have reduced computationally to the homogeneous equation ( 60 ) . 


With the procedures just outlined, the mesh points fall into eight 
categories, each of which requires a different treatment. These categories 
are listed as follows (a typical point in each category is indicated by a 
corresponding letter in sketch (d)): 

a. Subsonic points in q ,9 coordinates 

b. Regular supersonic points 

c. Supersonic points next to corner characteristic 

d. Supersonic points next to sonic line 

e. Points on sonic line 

f . Points near intersection of sonic line and corner characteristic 

g. Regular subsonic points in w,0 coordinates 

h. Points next to stagnation point 

In most categories the finite- difference approximation will be obtained 
by some variation of the usual series-expansion procedure. In this proce- 
dure the function F at each of a number of mesh points adjacent to the 
given point is expanded in terms of a truncated Taylor's series about the 
given point. These series can be regarded as simultaneous equations for 
the derivatives of F at the given point and solved accordingly. Follow- 
ing this approach, we obtain expressions for the derivatives in terms of 
the local dimensions of the mesh and the values of F at the points 


5 In the earlier work of references 3 and k, attempts to devise a 
numerical procedure in the supersonic region, using either rectangular or 
characteristic coordinates, were not successful. In these reports, which 
employed the Tricomi approximation, the difficulty was overcame by elimi- 
nating this region and substituting an integral relation as a boundary 
condition on the sonic line. For the present work it has been found pos- 
sible to set up a successful numerical procedure by having proper regard 
(which was not done previously) for the region of dependence of each point 
in the characteristics mesh. This was essential here, since an integral 
relation like that used before is not available for the exact equation. 
Even if it were, however, the present method would probably be preferred 
from the computational point of view. 
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considered. By substituting these results into the pertinent differential 
relations, the required finite-difference approximations are obtained. 

The number of terms retained in the Taylor’s series in any given case, as 
well as the number of adjacent points utilized, will depend upon the order 
of the derivatives to be approximated, the local arrangement of the mesh, 
and the degree of accuracy desired. The policy here will be to use the 
simplest expressions possible consistent with the analytical and geometri- 
cal situation at the point in question. The requisite accuracy will then 
be obtained by reduction of the mesh interval to whatever degree appears 
needed in various parts of the field. The details of the derivation will 
be given only when the procedure differs from that ordinarily encountered. 
(For detailed examples of the ordinary kind, see ref. 3*) 


Subsonic points in q ,9 coordinates .- The situation at points of 
this kind is shown in sketch (e) . For application next to the upper and 



lower boundaries, irregular mesh intervals are 
provided in the vertical direction. These are of 
length h and k relative to the basic interval 
A. By straightforward application of the series- 
expansion procedure to the approximation of the 
differential operator ( 5 ^a) , one obtains for the 
difference operator at 0 


l a (f) - A [; 


F2+F4 


kFi+hF 2 
1q hk(h+k) 




+ bcFo 


Sketch (e) 


(61) 


where the subscripts 0, 1, 2, etc., denote values of F at the correspond- 
ing points in the sketch. (The quantity b Q that appears in this and 
subsequent finite- difference expressions should not be confused with the 
bo that was introduced previously in the treatment of the free- stream 
singularity; cf. eqs. (21) through (ho) .) When the operator ( 6 l) is 
applied to i|r at a point 0 next to the upper or lower boundary, \lr 1 or 
is set equal to zero as required by the boundary conditions. For a point 
0 on the line q = q^, t 4 is replaced by T 4 /g 4 to conform with the 

variables used in the region q < q^ (see above). 


2 



0 


Sketch (f) 


Regular supersonic points .- The arrangement 
of points in this case is shown in sketch (f ) . 
Here the adjacent points 1, 2, 3 are chosen to 
lie in the region of dependence of the point 0 , 
which for the present flow is the curvilinear 
quadrant opening toward the corner characteristic 
and the sonic line.® Through application of 
the series-expansion procedure in the r,Z 


6 Note that the increment in r or I from one mesh line to the next 
is the same as A, the mesh interval on the sonic line (cf. eqs. ( 15 ) ) . 
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coordinates, the following difference operator is found "by approximation 
to the differential operator (5^b) : 


l a (f) 


[ 6 (r o -Z 0 )] 2/3 

A 2 



Fi - F3 + Fq) 


+ 


A 

6(r 0 - Z Q ) 



+ ^qFq 


(62) 


Supersonic points next to comer characteristic .- The situation at 
a point of this kind is shown in sketch (g). It is known from Guderley’s 
solution for flow over a convex comer (ref. l 6 ‘ 
that t will have a singularity in the second 
derivatives at the corner. This means that the 
series-expansion procedure is not permissible 
in this cane and some other approach must he 
devised. 

For values of A small enough to obtain 
acceptable results in the numerical solution, 
equation ( 59 ) will ordinarily have reduced to equation ( 60 ) by the time 
the boundary is reached. We need concern ourselves, therefore, only with 
an expression for L A as applied to i|r. According to Guderley * s singular 
solution, in the vicinity of the corner characteristic the variation of i|r 
along a line of constant Z (such as 03 or 12) is of the form 

* ~ ( r DE - r ) 4/S ( 6 3) 



where r DE is the value of r at the corner. Although Guderley* s solu- 
tion is for the Tricomi equation, the same result may be expected for the 
exact equation (7) . An interpolating function that satisfies condition ( 63 ) 
and may therefore be used as an approximation to iff in the vicinity of the 
point 0 is 



(l - 





Differentiation of this relation and substitution into expression (5^b) 
gives for the difference operator at the point 0 


%(♦) = 


1 6 ( r o ” l 0 ) ] 


2/3 


-jp (+1- + 


6 (r 0 - lo ) 


[C 1 ' *0 - +x]} +b °*o 

(64) 
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Sketch (h) 


Supersonic points next to sonic line .- 
A complication arises in this case (see 
sketch (h) ) from the fact that the deriva- 
tives of ^ in the r,l coordinates are 
singular on the sonic line (cf. eq. ( 17 )). 
As a consequence, a series expansion in 
the r ,2 coordinates, like that employed 
at regular supersonic points, cannot he 
•used here. This difficulty can he over- 
come hy using an expansion in tj,0 coor- 
dinates, -where the derivatives are regular, 
and approximating the differential operator 
in the form ( 5 ^a) . 

To carry out the approximation, the 
truncated Taylor's series for the func- 
tion F in the vicinity of the point 0 
is taken in the form 


F - F 0 + (n - n )F + (0 - + 


0) - n 0 )‘ 


Ft-t, + 


" 0 o) F TJ0 o + 


(0 - e o y 




In this connection we note that, for mesh points in the vicinity of 0, 

(tI-tIq) is of order A 2/3 whi le (0 - 0 O ) is of order A (see eqs. (l 6 )). 

This means that the cubic term in (q - t] q ) is of the same order in A as 

the squared term in (0 - 0 O ) and must he retained in the series . 7 By 

application of the series ( 65 ) successively at the points 1 through 'J , a 

set of simultaneous equations is obtained that can he solved for F__ and 

‘l'lo 

F qq q . Substitution of the resulting expressions into the differential 
operator ( 5 ^a) gives finally for the difference operator at the point 0 


Ia(f) 


1 

- /1 

(3A/4) 



+ F 4 ) + 0 (F 3 +F s ) + 2 tF s 




k 

+ ° + 3 


T 



+ b^o 


(66) 


7 This arises from the fact that, although the expansion is being 
carried out in the q ,0 variables, the mesh points are still arranged on 
lines of constant r and 2 . 
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•where 


Tv = 


2 2/S (2 - 2 W3 ) - 3 2/S (2 - 3 a/S ) 

g2/3^2/3 _ 2 2/3 ) 


cr = - 


2-3 


2/3 


22/3^22/ 3 _ l)(3 2/s - 2 a/3 ) 


T = 


2 - 2 ‘ 


2/3 


32 / 3 ^ 2/3 _ l)(32/3 _ 2 2/s ) 


When the point 0 is close to the corner - as at the point d* -in 
sketch (d) - a further complication arises since the points 5 , 6 , and 7 
of sketch (h) are no longer available. As will be seen later, however, 
analytical knowledge of the - nature of \|r where the sonic line approaches 
the corner characteristic can be used to obtain expressions for ^ and i|r 0 

at point 2. (Here only i|r is of interest, since at point d’ equa- 
tion ( 59 ) will have reduced again to equation ( 60 ) . ) By differentiation 
of the series ( 65 ) with respect to q and 9 and application of the result- 
ing expressions at 2 , simultaneous equations can again be obtained that 

can be solved for it and . In this way the following formula is 

yy o 

obtained for the necessary difference operator when point 0 is located 
at d»: 


I&(t) TvCti" t 2 ) + + T(i|r 4 -i|r 2 ) - 

(3AA) L 

i(x + 0+ T+|)(t 0 -t 2 ) +2x(f) + 

i (\ + a - 3t - l) Ai|r 0 J + b 0 T|r 0 


(67) 
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where now 


gg/3 + a i/3 _ a 
2 s ' 3 - 1 


CT 


, 1/3 


(a 273 - i) 


2 


T 



2 

1 



Points on sonic line .- The 
arrangement at a point of this type 
is shown in sketch (i) . The Taylor’s 
series for F as an expansion about 
point 0 is again taken in the form of 
equation ( 65 ) , where now t | 0 = 0 . 

By applying this series at points 1 
through 6 and solving the resulting 
equations for F__ , one can write the 
• *0 

difference operator, by approximation 
to equation ( 5 ^a) , as 


L ( F ) = _ i t(t 2 -l)F 2 - 2t(t 2 -4)F x - | (Fs+P 6 ) + 

A 2 t (t + 1 ) (t + 2 ) L 4 

3(F 4 + F 5 ) - ^ + 7t - t ^ F 0 J + boF 0 (68) 


where t is related to A by 


t = 



(69) 


points near intersection of sonic line and comer characteristic . - 
At the two points near the intersection of the sonic line and corner 
characteristic (points f x and f 2 in sketches (d) and (j) ) the difference 
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operators (68) anfl (66) cannot be used since the necessary adjacent points 
are not present. These two points can he eliminated from the problem, 
however, by use of an interpolating function in the region in question. 

An appropriate function is provided here by Guderley’s singular solution 
for flow over a convex corner (ref. l6) . Though determined for the Tricomi 
equation, this solution should provide a good approximation to the corre- 
sponding solution for the exact equation for small values of t). The 
details are as follows: 

According to Guderley’s solution, the variation of i|r in the vicinity 
of point D of sketch ( j) can be represented by 


^ - B|t|| fjj(£j2) 


( 70 ) 


where B is a constant and £ is now defined by 


9(9p - ef 

kt \ 3 


(71) 


The function fjj, given originally by Guderley in terms of hypergeometric 
series, can be represented in closed form by 


(£j2) - 



1 



1/3 



( 72 ) 


where now 


Qp-0 

A„= 



(73) 


To determine the constant B, the func- 
tion (70) is made to pass through the value 
of i|r at the mesh point 1 of sketch (j). 

We thus obtain finally for the interpolat- 
ing function in the region in question 


♦1 


S^-e ) 4 ' 3 


D 



|riffn:(5;2) (7*0 


Sketch ( j) 
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Thi s function is used 
wherever these values 
"boring points . After 
necessary expressions 
equation ( 67 ) . 


to eliminate the values of t|r for points f 1 and 
appear in the finite -difference equations for neigh- 
differentiation, equation (7*0 also provides the 
for the derivatives ijr and Tjr a required in 

^2 “2 



3 

Sketch (k) 


Regular subsonic points in w,9 coordinates . - 
A typical point, g, of this category is shown in 
sketch (k) . Here provision must he made for an 
irregular mesh interval in any direction from 0 as 
indicated in the sketch hy the constants e, f , h, 
s.nfl k. As previously explained, "the difference 
equation is obtained in this case by approximation 
to equation (l) rather than to equation (59) • ihe 
result, as obtained by the series-expansion proce- 
dure, can be written 


1 + 
f ( e 


where 


JL. y 2 + ¥ 4 +— — r (k^+Ma)- [ ~ e + f + (e = 

+ f) e(e+f) hk(h+k) L ef(e + f) hkj 


(75) 


P = 


1 - ^—4 w 0 2 

A 1 7+1 

^ w 0 2 

7 + 1 


q , (kt JL -l ua 
Q Va J w 0 2 . _ 


1 - 7 Wo‘ 

■7 + 1 


For application at points on the line T] = (see sketch (d)), Vg in 
equation ( 75 ) is replaced by g 2 i|r 2 as required by the variables used in 
the region tj > tj^. 

Points next to stagnation point .- The use of equation (75) a "k points 
next to the stagnation point - as, for example, at point h in sketch (d) 
is not practical. This is due to the fact that the final values of ¥ are 
obtained by taking the difference between two quantities which, according 
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to the boundary conditions at the stagnation point, are 
very nearly equal in this vicinity (see eqs. (3) 
through (5)). This difficulty can be overcome as follows: 

The situation in the region adjacent to the stag- 
nation point is shown in sketch (2). For small values of 
w, a solution of equation (l) that satisfies the boundary 
condition that IF = 0 on the surface of the plate (where 
Y . again represents either S 8, or T 8 ) is given by the 
series 


?(w,0) = Diw sin(e + a) + D 2 w 2 sin[2(0+ a)] + 

D 3 v={sin[3(0 + «)]-£ ♦ 0(w 4 ) (76) 



Sketch (l) 


Here a is the angle of attack and D x , D 2 , and D3 are constants. By 
writing equation (76) for any three points 1, 2, 3 on the vertical line 
w = 2A’ and neglecting terms of 0(w 4 ) , we obtain three simultaneous 
equations for D x , D 2 , and D3. These equations can be solved to obtain 




(77) 


and similarly for D 2 and D 3 . With these relations and equation (76) we 
can then express W at any point h on the line w = A* in the form 


\ ( 78 ) 

These expressions are used to replace Y4 in the difference equation (75) 
when that equation is written for mesh points on the line w = 2A' . In 
this manner points next to the stagnation point are eliminated from 
explicit consideration in the finite-difference scheme. For best accuracy 
the points 1, 2, 3 should be placed at approximately equal intervals 
between the upper and lower boundaries. If this is done, no difficulty 
is found in obtaining satisfactory results without a reduction in mesh 
size in the vicinity of the stagnation point. 

Distribution of mesh points .- The mesh used for an actual calculation 
is not uniform as in sketch (d) but has a different spacing in different 
parts of the field. The distribution of mesh interval for the present 
solution (a = 13°) is shown in figure 6. This particular distribution 
involves a total of ' 880 finite-difference equations. The same distribution 
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was used for ‘both ijr 0, and ijr®. As in references 3 and tlie transition 
between the various parts of the graded mesh was accomplished hy the use 
of overlapping fields in essentially the manner described in reference 18 . 


Solution of Finite-Difference Equations 


The solution of the finite-difference equations was obtained by an 
iterative process involving both relaxation techniques and step-by-step 
procedures. In the first stage of each iteration, the values of t for 
the column of mesh points immediately to the right of the sonic line (as 
obtained from an initial guess or from the previous iteration) are con- 
sidered fixed. On the basis of these values, a boundary- value problem 
is then solved for the points in the subsonic field and on the sonic line. 
This is done by the use of standard relaxation techniques (see, e.g., 
refs. 19, 20, and 21). In the second stage of the iteration the values 
of ijr obtained on the sonic line as a result of the first stage are used 
as the initial values for an initial-value problem in the supersonic field. 
This problem is solved in simple step-by-step fashion proceeding along 
the characteristics running from the sonic line to the limiting character- 
istic. This is done first for the characteristic adjacent to the corner 
nnrt then for succeeding characteristics in the interior of the flow. By 
tuphtih of this procedure one obtains a new set of values for at the 
points immediately to the right of the sonic line. Those values are now 
considered fixed again, and the entire process is repeated. This con- 
tinues until a consistent solution is obtained for the complete set of 
finite-difference equations . 

The numerical work for the present solution waB performed entirely 
on a desk calculator. In the initial stages of the iteration process the 
relaxation solution of the subsonic field was rarely carried to completion 
at any given stage. Over- and under- relaxation were, in fact, found use- 
ful to counteract the changes fed back by the subsequent work in the super- 
sonic region. Within the subsonic field, block relaxation was used exten- 
sively. Such devices were, in fact, essential to the solution of the 
problem in a practicable length of time. 

The results of the calculations are shown in figures 7 to 9 as contour 
maps of T in w ,0 coor dinat es. Numerical values are also given in 
table II at the end of the report. For the reasons expla in ed in part I, 
the calculations were made for a value of 7 of 1.405. The results for 
T 0, in figure 7 and table FT correspond to A a = 100; the results for T 8 
in figure 8 and table II to A s = 1,000. These values, which fix the 
strength of the singularities, were chosen to give a convenient level for 
the dependent variable throughout most of the field and to provide an 
absolute value of approximately 1 for the constant C in equation (^) . 

For both solutions the work was done to integer values of ijr in most of 
the field, and the residuals in the relaxation work were eliminated to 
within 12 (with due care that residuals in any given area were not pre- 
dominately of the same sign). Near the boundaries, where i|r is s mall . 
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the work was done with values of f to two decimal places, and the 
residuals were el iminat ed to within ±0.02. The final solution for 
given in figure 9 ar| d table II , was calculated from equation (4) with 
C = -1.007. Thi s figure was arrived at from the relation 


C 



which is obtained by substituting equation ( 76 ) into equation ( 5 ). The 
values of Di a and Dx s were computed from an equation of the form (77 ) , 
as obtained in the treatment of the mesh points next to the stagnation 
point. 


Transformation to Physical Plane 


The transformation from the holograph plane to the physical plane is 
accomplished by integration of equations (51) or (52) . Since the necessary 
integrations are chiefly along lines of constant w, 0 , r, or 2 , the equa- 
tions simplify considerably in application. The actual integration was 
carried out in the present work by mechanical means on the basis of plots 
of the appropriate integrand. The values of the derivatives of ¥ that 
appear in the various integrands were found by numerical differentiation 
of the values listed in table H. A detailed explanation of the procedure 
as applied to a similar problem has been given in reference 3 , and little 
more need be said here. It should be noted, however, that when the inte- 
gration of equations ( 52 ) is performed along a line of constant r or l, 
the exprepsion to be integrated reduces to the form f(w,0)df. In this 
case, therefore, no numerical differentiation of ¥ is required, which 
is advantageous with regard to the accuracy of the final result. 

With the distribution of w known on -the lower surface of the plate, 
the corresponding pressures can be calculated from the usual isentropic 
relations. This follows from the fact that the flow on the lower surface 
is irrotational (see Description of Flow Field). 


CHARACTERISTICS CONSTRUCTION OF FLOW 
OVER REGION OF SEPARATION 


As explained in part I, the calculation of the region of separation 
on the upper surface was accomplished on an electronic computer by means 
of a numerical method of characteristics. The method employed was based 
on the original Prandtl-Busemann procedure for two-dimensional irrotational 
flow (Bee, e.g., ref. 8). For the results shown in figure 1 the increment 
in stream angle between successive points in the characteristics net was 
taken as 1° throughout most of the field. In a small region adjacent to 
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the leading edge this Increment "was reduced to 0.2°. A check calculation 
using an increment of 2° outside the latter region gave results negligibly 
different from those obtained with 1°. With the 1° increment, the number 
of calculated points in the region of dependence of the separated stream- 
line was approximately 6600. To obtain the required initial data on the 
limiting wave (i.e., x/c and y/c as functions of r) it was necessary to 
interpolate between the values obtained from the mesh points used in the 
transonic calculation (see fig. 6). This was done graphically over most 
of the range. Near the leading edge, however, the interpolation was done 
analytically on the basis of an interpolating function compatible with 
the nature of the singularity in the vicinity of the corner characteristic. 
The values of the initial data are listed in table IH for anyone who may 
wish to calculate the flow over an airfoil with another type of upper 
surface . 


Ames Aeronautical laboratory 

National Advisory Committee 
Moffett Field, Calif., 


for Aeronautics 
May 7, 1956 
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APPENDIX 

NOTATION 
Primary Symbols 


a * 

A 

b 

c 

C 

C P 


critical speed (i.e.. speed at which speed of flow and speed 
of sound are equal) 

constant determining strength of singularity 

function of tj in differential equation for i|r (see eqs. (7) 
and (12)) 

constants in series expansion for b (see eq. (13)) 
chord of plate 

constant in superposition equation (see eqs . ( 4 ) and (5)) 

p — p 

pressure coefficient, — = — — 




constants in series solution for T in vicinity of stagna- 
tion point (see eqs. (76) and (77)) 


e,f,h,k 


^ j(C j n ) y 

f n(£; n ) 

F 


g 


length of irregular mesh intervals relative to that of basic 
interval 

functions appearing in singular solutions of Tricomi equation 
(see eqs. (19) and (70)) 


general symbol for variable in linear operator (may represent 
\Jr, M, or X ) 

function of w in transformation from I to t (see eqs. (6) 
and (11)) 


%> h i 

i 

I 

L*( ) 
%( ) 


undefined functions of t, 
index of summation (see eq. (22)) 
integral defined by equation ( 49 ) 
differential operator (see eq. ( 54 ) ) 
difference operator approximating L^( ) 
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m maximum value of Index i 

M q uant ity introduced in equation (56) as shorthand notation 

for quantities appearing in equations (36) and (39) 

n exponent of | tj | in singular solutions of Tricomi equation 

(see eqs. (19) and (70)) 

p static pressure 


Pi'* 1 ! 

<L 

r,l 


total pressure 

constants in solution of equation for (see eq. (30)) 
dynamic pressure 

characteristic coordinates (see eqs. (15) and (l6)) 

any length in flow field; also function of tj and 0 defined 
by equation (4-2) 


(sign tj) 3 1 quantities with absolute value 1 and sign of n or 0 
(sign 0) J 

t index of summation (see eq. (26)) 

T( ) Tricomi operator (see eq. (l8)) 


w speed of flow made dimensionless through division by critical 

speed 

x,y Cartesian coordinates (x in direction of free stream, y 

vertical, origin at leading edge of plate) 


x» 


distance along plate measured aft from leading edge 


a angle of attack 

P function of w (see eq. (53)) 

7 ratio of specific heats (1.405 in numerical work) 



A 

A 1 


dimensionless distance in free stream between given stream- 
line and stagnation streamline 

basic mesh interval in tj and 0 

basic mesh interval in w 



function of rj and 0 defined lay equation (20) 

function of w in transformation from I to f (see eq. (10)) 

in cl in ation of flow measured counterclockwise from free- 
stream direction 

density of fluid 

total density of fluid 

function appearing in solution for \|r (see eq. (22)) 

•transformed stream function (see eq. (6)) 

solution representing free-stream singularity for Tricomi 
equation (see eq. (19)) 

singular correction functions (see eq. (22)) 

function appearing in solution for (see eqs. (27), (28), 
and (30)) 

stream function 

Subscripts 

conditions in free stream 
conditions on separated streamline 

values corresponding to maximum Bpeed (i.e. , to speed 
attained by expansion to a vacuum) 

values of rj and w at change from tj to w coor dinat e in 
finite-difference mesh 

value at prescribed mesh point 

Superscripts 

quantities pertaining to solution with antisymmetric 
singularity 

quantities pertaining to solution with symmetric s ingular ity 

first and second derivatives of function with respect to its 
argument 



MCA TN 3723 


48 


REFERENCES 


1. Guderley, Gottfried: The Flow Over a Flat Elate With a Small Angle 

of Attack at Mach Number 1. Jour. Aero. Sci., vol. 21, no. 4, 

Apr. 1954, pp. 261-274. 

2. Fr ank! , F.: Investigation of the Theory of Wings of Infinite Span. 

Moving at the Speed of Sound. Doklady Akad. Nauk SSSR, vol. 57, 
no. 7, 1947, pp. 661-664. 

3. Vincenti, Walter G. , and Wagoner, Cleo B.: Transonic Flow Past a 

Wedge Profile With Detached Bow Wave. MCA Rep. IO95, 1952. 
(Supersedes MCA TN*s 2339 and 2588) 

4. Vincenti, Walter G., and Wagoner, Cleo B.: Theoretical Study of the 

Transonic Lif t of a Double-Wedge Profile With Detached Bow Wave. 

MCA Rep. 1180, 1954. (Supersedes MCA TN 2832) 

5. Mepmann, H. W., and Bryson, A. E., Jr.: Transonic Flow Past Wedge 

Sections. Jour. Aero. Sci., vol. 17, no. 12, Dec. 1950 , PP. 745-755 
(Also published as C.I.T., Guggenheim Aero. Lab. Pub. 267) 

6. Guderley, Gottfried: Two-Dimensional Flow Patterns with a Free- 

Stream Mach Number dose to One. USAF, WADC, AF Tech. Rep. No. 6343 
May 195L. 

7. Guderley, G., and Toshihara, H. : The Flow over a Wedge Profile at 

Mach Number 1. Jour. Aero. Sci., vol. 17, no. 11, Nov. 1950, 

PP. 723-735. 

8. Guderley, G.: Monograph on the Theory of Characteristics. Transla- 

tion No. RAT-4, Douglas Aircraft Co., Inc., Santa Monica, Calif., 
Sept. 1, 1947. 

9. Durand, W. F.: Fluid Mechanics, Part I. Div. B of Aerody nam ic 

Theory, vol. I, W. F. Durand, ed., Julius Springer (Berlin), 1934, 
PP. 159, 184-185. 

10. Courant, R., and Friedrichs, K. 0.: Supersonic Flow and Shock Waves. 

Inters cience Pub., 3hc., 1948. 

11. Bickley, W. G. : Some Exact Solutions of the Equations of Steady 

Homentropic Flow of an InviBcid Gas. Chap. V of Modern Developments 
In Fluid Dynamics; High Speed Flow, vol. I, L. Howarth, ed., 

Oxford Univ. Press, 1953. 

Guderley, Gottfried: Transonic Simplifications of the Eodo graph 

Equation. WADC Tech. Rep. 53-183, Wright-Patterson Air Force 
Base, June 1953. 


12 . 





NACA TN 3723 49 


13. Wood, George P.: Experiments on Transonic Flow Around Wedges. 

MCA TN 2829, Not. 1952. 

14. Sauer, Robert: Introduction to Theoretical Gas Dynamics. Trans. 

by Freeman K. Hill and Ralph A. Alpher. J. W. Edwards, Ann Arbor, 

I9V7. 

15. Guderley, Gottfried, and Yoshihara, Hideo: Two-Dimensional. Unsym- 

metric Flow Patterns at Mach Number 1. Jour. Aero. Sci., vol. 20, 
no. 11, Nov. 1953, pp. 757-768. 

16. Guderley, K. Gottfried: Singularities at the Sonic Velocity. Tech. 

Rep. F-TR-II 7 I-ND, Air Materiel Command. Wright-Patterson Air 
Force Base, Tech. Intelligence, June 1948 . 

17. Woods, L. C.: The Relaxation Treatment of Singular Points in 

Poisson’s Equation. Quar. Jour. Mech. and Appl. Math., vol. VI, 
pt. 2, June 1953, PP. I 63 -I 85 . 

18. Tasny-Tschiassny, L. : The Triangulation of a Two-Dimensional 

Continuum for the Purpose of the Approximate Solution of Second- 
Order Partial Differential Equations . Jour. App. Fhys., vol. 20, 
no. 5 , May 1949, PP. 419-424. 

19 . Fox, L.: A Short Account of Relaxation Methods. Quart. Jour. Mech. 

and Appl. Math., vol. I, pt. 3, Sept. 1948, pp. 253-280. 

20. Shaw, F. S.: An Introduction to Relaxation Methods. Dover Pub., 

Inc., 1953. 

21. Allen. D. N. de G.: Relaxation Methods. McGraw-Hill Book Co., Inc., 

1954. 



50 


MCA TN 3723 


TABLE I.- VALUES OF v, tj, g, AND b (7 = 1.405) 


V 

n 

g 

-b 

0 

-00 

00 

0 

.109 

-2.037 

1.736 


.164 

-I.745 

I.671 


.218 

-1.524 

1.615 


.272 

-1.343 

1.566 

— 

.326 

-1.190 

1.520 

— - 

•379 

-1.056 

1.477 

— 

.432 

-•935 

1.436 

.261 

.484 

-.826 

1*395 

.300 

•539 

-.725 

1.356 

.343 

• 585 

-.631 

1.317 

.388 

•635 

-.545 

1.280 

.442 

.684 

-.463 

1.243 

.482 

•732 

-.386 

1.206 

.540 

•779 

-313 

1.170 

•593 

.825 

-.245 

1.135 

.652 

.871 

-.179 

1.100 

.712 

.915 

-.116 

1.065 

• 777 

• 958 

-.057 

1.035 

.850 

1.000 

0 

1.000 

.918 

1.067 

.088 

.947 

l.o4i 

1.107 

•l4o 

.917 

1.125 

i.l4i 

.183 

.890 

1.200 

1.200 

.258 

.842 

1.339 

1.252 

-323 

.800 

1.475 

1.323 

.409 

.744 

1.685 

1.428 

• 536 

.660 

2.055 

1.522 

.650 

• 585 

2.531 

1.609 

.754 

.515 

3-045 

1.690 

.851 

.452 

3-641 

1-765 

.943 

.394 

4.3IJ8 

1.836 

1.031 

.340 

5-193 

1.902 

1.115 

.291 

6.378 

1.964 

1.197 

.246 

7-461 

2.022 

1.275 

.204 

8.986 

2.076 

1-351 

.170 

10.894 

2.127 

1.425 

.138 

13.301 

i 2.173 

1.497 

.109 

l6.4o6 

2.215 

1.568 

.084 

20.468 

2.254 

1.637 

.065 

25.904 

2.289 

1.704 

.0 48 

33.394 

2.389 

1-950 

.oil 

- — 

2.437 

2.252 

0 

00 
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TABLE H.- T 8 -, 5 s , AND 7 AS FUNCTIONS OF v AND 6 (A a = 100, A s = 1000, 

7 = 1A05) 
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IABIE H.- 5 s , AMD f AS FUNCTIONS OF v AND 9 (A a = 100, A s = 1000, 

7 = 1.405) - Continued 


e 

s 4 

I 8 

T 

V 

B 

s* 

I s 

T 

w 

B 

5 s 

J 8 

S 

0.10 

3,5&) 

3,561 

-8 

0.604 

0.30 

1,792 

1,596 

184 

0.450 

0.70 

a? 

3® 

30 

.05 

2,740 

3,899 

-1 A® 


•25 

1,948 

1,794 

l4i 


.60 

704 

670 

29 

0 

1,165 

3,883 

-2,746 


.20 

2,044 

1,992 

38 


-5° 

793 

767 

20 

-.05 

-375 

3,423 



J5 


2,169 

-158 


.40 

867 

862 

-2 

-.10 

-1,120 

2,599 

-3,727 


JO 

1,844 

2,289 

-462 


.30 

901 

943 

-50 

-.15 

-997 

1,606 

-2,614 


.05 

1,478 

2,308 

-848 


.20 

867 

978 

-118 

-.20 


564 

-967 


0 

990 

2,1® 

-1,211 


.10 

740 

928 

-195 

2.80 

2.70 

6.49 

12.31 

6.26 

11.88 



-.05 

-JO 

-.15 

510 

164 

1 

m 

-l,4o4 

-1,312 

.405 

1 O 

e 

535 

296 


-235 

-188 

2.60 

18.41 

17.75 

.672 


-.20 

-16 

328 

-346 

2.® 

10.72 

10.32 

•331 

2.50 

24.87 

23.94 

.924 


.45 





2.® 

29.80 

28.® 

.984 

•90 

423 

38l 

40 

.576 

1,210 

1,071 

131 


2.40 

50.44 

48.51 

1-71 

.80 

528 

467 

57 


.40 

E 

1A99 

135 


2.20 

73-8 

70.9 

2.58 

.70 

674 

5® 

® 


.35 

E'El 

1,342 

130 


2.00 

101.2 

96.9 

3-® 

.60 

878 

745 

129 


.30 

1,609 

ml 1 

106 


1.80 

134.6 

128.6 

5.15 

.50 

1,172 

970 

194 


.25 

1,708 

If ^.7l 

51 


1.® 

176a 

167.8 

7-1 

>5 

1,363 

1,119 

237 


.20 

1,750 

1,790 

-53 


i.4o 

229.2 

216.0 

9-7 

.40 

1,593 

1,297 

287 


-15 

!,704 

1,908 

-ZL7 


1.20 

£97-9 


12.8 

•35 

1,8® 

1,512 

341 


JO 

1,544 

1,969 

-44o 


1.® 

386.5 


16.2 

.30 

2,177 

1,768 

397 


.05 

1,264 

1,946 

-697 


•90 

*39 


16 

.53 

2,515 

2,071 

428 


0 

913 

1,819 

-919 


.® 

497 

477 

16 

.20 

Kl J 

2,412 

389 


-05 

5® 

1,571 

-1,009 


.70 

560 

541 

15 

-15 


2,768 

1® 


-.10 


1,208 

-933 


.60 

623 

608 

10 

.10 

2,801 

3,087 

-309 


-J5 

117 

757 

-644 


•50 

678 

674 

-1 

•05 

2,152 

3,270 

-1,142 


-.20 

32 

2® 

-236 


.40 

714 

730 

-21 

0 

1,121 

3,193 

-2,094 



8.34 




•30 

7?9 

767 

-54 

-.03 

U6 

2,810 

-2,714 

-550 

2.® 

8.03 

.292 


.20 

®5 

‘763 

-93 

-.10 

-460 

2,160 

-2,635 


2.70 

15.75 

15.16 

-555 


.10 

573 

697 

-129 

-.15 

-517 

1,353 

-1,880 


2.® 

23-39 

22-53 

.835 


0 

24l 

557 

-l4o 

-;P0 

-221 

478 

-702 


2-50 

•90 

31-39 

467 

30.28 

429 

1.138 


-.10 

344 

-107 

.45 

1,328 

1,116 

204 


.80 

5® 

514 

47 

.362 

.90 

412 

4® 

9 

.40 

1,523 

1,280 

236 


•70 

689 

621 

63 


.80 

458 

446 

8 

-35 

1,749 

l,4g 

2 66 


.60 

843 

757 

81 


.70 

504 

495 

5 

-30 

1,98 9 


28a 


.50 

E 53 

931 

94 


.60 

546 

542 

1 

-25 

2,223 

1*942 

267 


.40 


1,141 

?4 


.50 

576 

5® 

-13 

.20 

2,403 

2,207 

181 


•30 

1,434 

1,377 

47 


•4o 

591 

613 

-27 

.15 

2,445 

2,462 

-34 


.20 

1,504 

1,594 



.30 

577 

622 

-49 

.10 

2,248 

2,6® 

-436 


JO 

1,307 

1,6® 

-394 


.20 

534 

598 

-69 

-05 

1,760 

2,748 

-1,007 


0 

831 

1,517 

-697 


JO 

451 

530 

-83 

0 

1,062 

2,639 



-JO 

336 

999 

-670 


0 

333 

415 

-® 

-.05 

381 

2,311 


.498 






-.10 

194 

254 

-62 

-.10 

-63 

1,7® 

-1,860 

2.® 

9.21 

8.® 

.318 


2.® 




-.15 

-198 

1,125 

-1,330 


2.® 

25.71 

24.73 

.900 

.323 

11.75 

11-33 

.355 

-.20 

-96 

399 

-498 


2.40 

43. 81 

42.08 

1.58 


2.® 

32.57 

31.40 

•990 



2.20 

64.6 

62.0 

2.43 


2.40 

54.8L 

52.81 

1.® 

2.80 

7.42 

7-15 

.263 


2.00 

®.6 

®.6 

3.57 


2.20 

79-45 

76.48 

2.50 

2.70 

l4.o4 

13-53 

.501 


1.® 

120.7 

114.9 

5.16 


2.00 

107.63 

103.50 

3.47 

2.60 

20.90 

20.14 

.756 


1.® 

161.1 

152.5 

7.5 


1.80 

i4o.to 

135.16 

4.62 

2.50 

28.15 

27.ll 

1.036 


1.4o 

215.3 

202.7 

n.i 


1.60 

i®.15 

172.95 

5.95 

2.40 

35.90 

34.50 

1.349 


1.20 

290.7 

271.9 

16.8 


1.40 

227.57 

218.® 

7.38 

2.20 

53.4 

51.2 

2.11 


1.® 

398.4 

370.4 

25.4 


1.® 

284.13 

273.® 

8.51 

2.00 

74.7 

71.2 

3.14 


•90 

469 

436 

30 


1.® 

349.27 

338.43 

8.40 

1.80 

■Exa 

96.5 

4.64 


• ® 

556 

513 

39 


.® 

418.06 

409.99 

5.12 

1.60 

137.5 

129-7 

6-9 


•70 

6® 

609 

47 


.® 

476.29 

477-25 

-4.4o 

1.4o 

187.4 

175.2 

10.9 


-® 

7® 

723 

50. 


.40 

493-70 

513.63 

-23.63 

1.20 

260.1 

240.6 

17.8 


•50 

913 

856 

51 


.20 

430.05 

475.01 

-48,39 

1.00 

372.2 

339-7 

30.1 


.40 

1,043 

1,001 

34 


0 

264.98 

315.38 

-52.67 

.90 

450 

409 

38 


•30 

1,134 

W*9 

-24 






.80 


497 

54 


.20 

1,127 

p.j 

-133 

.200 

2.® 

11.28 

10.94 

.262 

•70 

6p* 

613 

76 


.10 

968 

E'ja 

-286 


2.® 

31.09 

30.16 

.720 

.60 

879 

766 

108 


0 

672 

1,068 

-4o4 


2.40 

51.57 

50.04 

1.18 

•50 

1,126 

970 

150 


-.10 

343 

686 

-349 


2.20 

73.12 

70.97 

1.® 

•45 

1,273 

1,098 

168 

.450 


4® 

433 



2.00 

96.06 

93-30 

2.10 

•4o 

1,437 

1,245 

184 

•90 

24 


1.® 

120.59 

117.27 

2.48 

•35 

1,614 

i,4n 

193 


• ® 

532 

501 

28 


1.® 

146.59 

142.® 

2.71 
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TABLE II.- S 8 , AMD Y AS FUNCTIONS CF w AMD 0 (A a = 100, A s = 1000, 

y = 1.405) - Continued 


V 

e 

2® 

7“ 

T 

V 

8 


3 s 

I 

V 

6 

S* 

2® 

2 

0*200 

fw 

173-48 

169.61 

2.65 

1.21B 

2.80 

..o496 

.0483 

.00178 


•275 


266 

— 7 T 


1.20 

199.87 

196.35 

2.11 


2.70 

[1.024 

•997 

•0352 


.225 

545 

319 

■ 106 


1.00 

223.12 

220.72 

.811 


2.60 

2'4o6 

2.339 

.0843 


.175 

389 

153 


.80 

238.92 

238-72 

-1.52 


2.50 

4.077 

3.947 

.1454 





.ooi4o 


.60 

241.05 

2W-.25 

-4.96 


2.40 

5.961 

5.757 

• 2173 

1-357 

2.70 

.0391 

.0380 


.40 

222.09 

229.31 

-8.87 


2.30 

8.005 

7-724 

• 2997 


2.60 

.788 

.767 

.0270 


.20 

176.14 

186.15 

-U.35 


2.20 

10.26 

9.88 

• 392 


2.5O 

1.851 

1.797 

.0648 


0 

103.42 

111.87 

-9.255 


2.10 

12.64 

12.16 

• 502 


2.4o 

3-127 

3.028 

.1116 


.425 


447 



2.00 

15-25 

14.59 

.623 


2.30 

4.558 

4.403 

.1663 

1.083 

578 

128 


1.90 

18.18 

17.34 

.767 


2.20 

6.096 

5-883 

.2286 


•375 

715 

535 

176 


1.80 

21.41 

20.42 

.922 


2.10 

7.78 

7.50 

.298 


• 335 

908 

650 

254 


1.70 

24.95 

23.67 

1.139 


2.00 

9.54 

9.18 

.381 


■275 

1,182 

807 

369 


1.60 

28.96 

27-33 

1.38 


1.90 

11.45 

10.95 

.470 


.223 

1,610 

1,024 

579 


1.50 

33.72 

31.72 

1.71 


1.80 

13.59 

12.93 

.574 


•175 

2,339 

1,352 

978 


1.40 

39.32 

36.8O 

2.14 


1.70 

15.93 

15.16 

.698 


•125 

3,759 

1,878 

1,867 


1.30 

45.81 

42.80 

2.70 


1.60 

18.47 

17.49 

.845 


• 075 

7,229 

2,882 

4,326 


1.20 

53.81 

49.92 

3-48 


I.50 

21.35 

20.11 

1.03 


.025 

21,717 

5,456 

16,221 


1.10 

63.71 

58.82 

4.4? 


1.S0 

24.75 

23.28 

1.26 


-.025 

-20,119 

4,83?. 

24,993 


1.00 

76.57 

70.02 

6.0S 


1.30 

28.76 

26.92 

1-59 


-075 

-5,513 

2,o44 

-7,572 


•90 

92.18 

83.34 

8.25 


1.20 

33.44 

31.22 

2.00 


-125 

-1,915 

853 

-2,775 


.80 

n4.6 

103.0 

10.8 


1.10 

39.22 

36.35 

2.57 


-.175 

-4l4 

184 

-598 


.70 

144.5 

128.3 

15.2 


1.00 

46.4o 

42.79 

3.30 







.60 

188.0 

163.2 

23.6 


.90 

55-79 

50.94 

4.48 

1.138 

2.85 

•0579 

.0564 

.00208 


.50 

253.8 

211,8 

4o.5 


.80 

67-23 

60.66 

6.13 


2.75 

1.215 

I.I83 

-04lB 


.So 

357 

285 

70 


.70 

83.8 

75.2 

8.0 


2.65 

2.847 

2.765 

.0997 


•35 

433 

335 

95 


.60 

106.1 

9^.0 

H.5 


2.55 

4.809 

4.65? 

.1716 


■30 

534 

4oi 

131 


•50 

138.4 

119.9 

17.7 


2.4? 

7-026 

6.784 

.2560 


•25 

679 

487 

189 


.4o 

187-9 

156.2 



2.35 

9.433 

9.102 

• 3533 


.20 

883 

6o4 

276 

435 


.30 

266 

210 


2.25 

12.10 

U.65 

• 462 


•15 

1,206 

766 


.25 

322 

248 

73 


2.15 

14.93 

14-3? 

.593 


.10 

1,750 

1,008 

734 


.20 

4oi 

298 

101 


2.05 

1B.03 

17-24 

.736 


-.10 

-306 

136 

4% 







1-95 

21.53 

20.55 

.907 



4oo 



1.387 

.225 

301 

231 

68 


1.85 

25-42 

24.25 

1.113 

1.257 

•325 

309 

88 

1.417 






1-75 

29.67 

28.19 

1.355 


•275 

495 

371 

122 

2.65 

• 0352 

.0342 

.00126 


I.65 

34.52 

32.59 

1.65 


.225 

629 

451 

175 


2.55 

• 703 

.684 

.0242 


1.55 

40.26 

37.91 

2.o4 


■175 

820 

560 

256 


2.45 

1.648 

1.601 

.0577 


1.45 

47.03 

44.02 

2.56 


•125 

1,119 

709 

4o6 


2.33 

2.778 

2.690 

.0992 


1.35 

54.86 

51.24 

3.22 






2.25 

4.037 

3-901 

•i4t4 


1.25 

64.51 

59-86 

4.i4 

1.292 

2-75 

.0436 

.0424 

.00156 


2.15 

5.385 

5.198 

.2020 


1.15 

76.39 

70.56 

5.30 


2.65 

.889 

.866 

.0305 


2.05 

6.8 6 

6.6O 

.263 


1.05 

91.84 

84.02 

7.21 


2.55 

2.091 

2.Q31 

.0731 


i.95 

8.38 

8.06 

•335 


•95 

110.55 

100.02 

9.80 


2.45 

3.538 

3-426 

.1263 


1.85 

10.03 

9.57 

.411 


•85 

137-3 

123-7 

12.7 


2.35 

5-169 

4.993 

.1884 


1.75 

U.87 

H.29 

.500 


• 75 

173.1 

153.7 

iB.4 


2.25 

6.93 

6.688 

.2597 


1.65 

13.87 

13.18 

.608 


•65 

224.7 

195.3 

£8.0 


2.15 

8.87 

8.54 

.339 


1.55 

16.04 

15.19 

•736 


-55 

303.0 

253.1 

48.1 


2.05 

10.90 

10.49 

.434 


1.45 

18.48 

17-42 

.890 


• 49 

425 

339 

83 


1.95 

13-12 

12.55 

• 538 


1.35 

21.38 

20.12 

1.10 


.40 

516 

399 

u4 


I.85 

15-61 

14.87 

•639 


1.25 

24.83 

23-23 

1-38 


• 35 

636 

478 

155 


1.75 

18.34 

17.46 

.8o4 


1.15 

28.82 

26.89 

1.73 


.30 

809 

580 

2® 


1.65 

21.32 

20.22 

■ 974 


1.05 

33-78 

31.29 

2.23 


•25 

1,050 

719 

326 


1-55 

24.70 

23.28 

1.18 


• 95 

39-96 

36.83 

2.88 


.20 

1,433 

911 

516 


1.4? 

28.70 

26.99 

1.46 


.85 

48.08 

43.83 

3-91 


•15 

2,078 

1,198 

871 


1.35 

33.41 

31.28 

I.83 


• 75 

57.98 

52.26 

5-35 


.10 

3,319 

1,667 

1,640 


1.25 

38.87 

36.30 

2.31 


•<5 

72.4 

65.O 

6.9 


•05 

6,370 

2,550 

3,801 


1.15 

45.64 

42.32 

2.97 


• 55 

91.8 

81.2 

10.0 


-.05 

-4,824 

1,794 

-6,630 


1.05 

54.01 

49.84 

3.80 


M 

120.2 

103-9 

15-6 


-.10 

-1,666 

742 

-2,4l4 


•95 . 

64.91 

59.31 

5.17 


•35 

163.4 

135.6 

26.8 


-l 5 

-364 

161 

-526 


.85 

78.17 

70.62 

7.04 


•25 

232 

183 

48 



469 




• 75 

97-2 

87.4 

9.2 

1.470 





1.179 

•375 

364 

103 


.65 

123.0 

109.0 

13-1 

2.60 

.0319 

•0310 

.00115 


•323 

581 

43 6 

142 


.55 

160.0 

138.9 

20.2 


2.50 

.632 

.615 

.0217 


•275 

736 

529 

204 


.45 

216.6 

180.3 

34.8 


2.40 

1.477 

1.435 

.0517 


.225 

958 

65 6 

298 


• 35 

305 

243 

60 


2.30 

2.481 

2.403 

.0887 


•175 

1,307 

830 

471 


.30 


286 

82 


2.20 

3.596 

3-475 

.1313 


.125 

1,895 

1,094 

794 


•25 

|70 

584 

343 

113 


2.10 

4.781 

S.615 

.1795 


•075 

3,017 

1,520 

1,486 


.20 

418 

163 


2.00 

6.07 

5.84 

•233 


-075 

-1,505 

675 

-2,185 


.15 

761 

518 

239 


1.90 

7.39 

7-10 

.295 


-Mg 

-331 

147 

-478 







1.80 

8.82 

8.42 

.361 
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TABLE II.- ¥ a , f 6 , AM) H AS FUNCTIONS CF w AND 0 (A a = 100, A s = 1000, 

7 = 1.405) - Concluded 


V 

0 

D 

I s 

¥ 

V 

0 


7 s 

¥ 

1.914 

1.00 

2.92 

2.75 

■ BEi 

2.077 

MESM 

.0070 

.0068 

.00025 


.90 

3-23 

3.03 

.171 


1.75 

.116 

•113 

.oo 4 i 


.80 

3.65 

3 . 4-1 

.211 


1.65 

.249 

.242 

.0088 







1.55 

.382 

.371 

.0138 

1.950 

2.05 

.0110 

.0107 

.00039 


1.45 

•504 

.487 

.0187 


1.95 

• 193 

.188 

.0067 


1.35 

.605 

.584 

.0232 


1.85 

.428 

. 4-16 

.0151 


1.25 

.69 

•67 

.0270 


1-75 

.680 

.659 

.0245 


1.15 

• 76 

• 73 

.0318 


1.65 

• 931 

.900 

.0343 


1.05 

.81 

•77 

• 0359 


| 5 >/| 

1.166 

1 . 124 - 

.0443 








1-39 

1 . 34 - 

.0540 

2.106 

1.80 

.0061 

.0059 



1-35 

1-59 

1-53 

.0650 


1.70 

.100 

.097 

.0035 


1.25 

I.78 

1.70 

•0758 


1.60 

.211 

.206 

•P075 


1.15 

1.99 

1.88 

.0878 


1.50 

.321 

■ 311 

.0116 


1.05 

2.19 

2.08 

.103 


i. 4 o 

. 4 i 8 

. 4 o 4 

.0155 


• 95 

2.40 

2.27 

.121 

- 

1.30 

.496 

.478 

.0191 


.85 

2.65 

2.48 

.143 


1.20 

•56 

•54 

.0220 







1.10 

.60 

.58 

.0257 

1.983 

2.00 

.0099 

.0096 

.00036 







1.90 

.172 

.168 

.0060 

2.133 

1.75 

.0054 

.0052 

.00019 


1.80 

.378 

.368 

.0134 


1.65 

.085 

.083 

.0030 


1.70 

.597 

• 579 

.0215 


1.55 

.178 

.174 

.0063 


1.60 

.811 

• 784 

.0299 


1.45 

.268 

.260 

.0097 


1.50 

1.007 

.971 

.0384 


1.35 

.344 

-333 

.0128 


1 . 4 o 

1.19 

1 . 14 - 

.0463 


1.25 

.401 

.387 

.0155 


1.30 

1.35 

1.29 

.0552 


1.15 

.444 

•43 

.0178 


1.20 

1.50 

1 . 4-3 

.0641 







1.10 

1.65 

1.57 

•0737 

2.161 

1.70 

.0046 

.0045 

.00017 


1.00 

1.81 

1.72 

.0861 


1.60 

.072 

.070 

.0025 


.90 

1.97 

1.86 

.101 


1.50 

.149 

.145 

.0053 







1 . 4 o 

.214 

.214 

.0080 

, 






1.30 

.278 

.269 

.oio 4 

2.016 

1*95 

.0088 

.0086 

.00032 


1.20 

.318 

.307 

.0124 


1.85 

.151 

.147 

.0053 







1-75 

.330 

.321 

.0117 

2.187 

1.65 

• 0040 

.0038 

.00014 


1.65 

.516 

.501 

.0187 


1-55 

.060 

.059 

.0021 


1-55 

•695 

.672 

.0257 


1.45 

.124 

.120 

.oo 44 


1A5 

.855 

.824 

.0326 


1-35 

.179 

.174 

.0065 


1.35 

1.00 

.96 

.0388 


1.25 

.223 

.215 

.0083 


1.25 

1.12 

1.08 

.0462 







1.15 

1.23 

1.17 

•0532 

2.211 

1.60 

.0034 

.0034 

.00012 


1.05 

!-35 

1.28 

.0608 


1.50 

.051 

.049 

.0018 


• 95 

1.46 

1.39 

.0707 


i. 4 o 

.102 

.099 

.0036 







1.30 

.145 

•i 4 i 

.0053 

2.046 

1.90 

• 0079 

.0077 

.00028 







1.80 

•133 

.130 

1.70 

.288 

.281 

1.60 

.447 

.434 

1.50 

•596 

•576 

1.40 

.726 

•699 

1.30 

.84 

.81 

1.20 

• 93 

.89 

1.10 

1.01 

• 96 

1.00 

1.09 

i.o 4 


.0047 

2.234 

i .55 

.0102 


1.45 

.0162 


1-35 

.0221 

.0278 

.0327 

2 i 257 

1.50 

b 


i. 4 o 

.0500 

2.280 



0029 

.0028 

.00011 

.o 4 i 

.o 4 o 

.0015 

.081 

.079 

.0029 

0024 

.0023 

.00009 

.033 

.032 

.0012 

0020 

.0019 

.00007 
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TABLE III.- x/c AMD y/c AS FUNCTIONS CO? r ON LIMITING CHARACTERISTIC . 


r 

(x/c) 

(y/c) 

2.9 

-0.087x10"® 

1. 667XL0"® 

2.8 

-. 665x10“® 

21.40X10“® 

2.7 

-.552x10"® 

44.60x10"® 

2.6 

•394x10"® 

67.20x10“® 

2.5 

2.127X10"® 

88.57x10“® 

2.4 

4.648X10“® 

110.15x10“® 

2.3 

7-911x10"® 

131.6x10-® 

2.2 

12. 53x10"® 

155.8x10-® 

2.1 

18.02x10“® 

179.9x10"® 

2.0 

24.97x10"® 

206.3x10"® 

1.9 

34.25x10"® 

237.3x10“® 

1.8 

46. 76x10“® 

273.9x10"® 

1-7 

62.54x10"® 

315-4x10"® 

1.6 

85.45x10“® 

371.4x10"® 

1.5 

118. 18X10“® 

443.9x10"® 

1.4 

I61.65XLO"® 

533.7XLO"® 

1.3 

223.73x10"® 

653.6x10"® 

1.2 

301.87XL0"® 

791-3x10"® 

1.1 

433.60x10"® 

1008.1x10"® 

1.0 

6l6 . 2XL0" a 

1298.4X10"® 

• 9 

8.330XL0" 4 

16.359XLO" 4 

.8 

12.05XL0" 4 

21.933XL0" 4 

•7 

1S.52XL0" 4 

31.01XL0" 4 

-6 

31.38XL0" 4 

49.Q2X10" 4 

•5 

55-68XL0" 4 

81.91X10" 4 

-45 

74-34XL0" 4 

107. 0X10" 4 

.4 

102.59X10" 4 

145.4X10" 4 

-35 

144.01XL0" 4 

202.6XL0" 4 

-3 

208.QXL0" 4 

292.4XL0" 4 

-25 

324.2XL0" 4 

457.1X10" 4 

.2 

529-4XL0" 4 

766.0X10" 4 

— . 

957.7XL0" 4 

1444. 7XL0" 4 
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(b) Vicinity of leading edge 
Figure 2.- Concluded. 
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Figure 5 .- "b . as a function of 
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